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On a Two-Dimensional Free-Boundary Problem 


By Ravpu J. BEAN, Ropert P. GrirBert? ) and Hatcoms KENDALL, 
Pittsburgh, Pa., U.S.A.2) 


1. Introduction 


In a petroleum reservoir the gas-oil-water interfaces frequently deviate 
from the horizontal. The factors causing this phenomenon are of two basic 
types, static factors such as variation in formation properties, and dynamic 
factors such as hydraulic circulation of the oil and the underlying water in the 
stratum [1, 2]*). Occasionally circulation is sufficient to even cause the oil-cap 
to migrate to a new position. In fact it has been suggested that one might 


_ artificially set in motion the water underlying an oil-cap in order to move the 


bho AS 


ee ee SN Ck 


oil to a position more favorable for production [3]. Consequently, it is of 
importance to the petroleum industry to be able to determine the position and 
shape of an oil-water interface under dynamic conditions. This is the problem 
investigated in this work. 

It would be well to list here the symbols used in the following work: 


Variables 


x,y distance coordinates 


t time 

Vo oil velocity 

v water velocity 

dbo oil scalar potential 

db, water scalar potential 

Py oil hydrodynamic pressure (dynes cm~*) 
P, water pressure 

Ve the stream function 


w=®+iYW the complex potential 


re oe zi Ve Ge Vy the complex ‘velocity’ 
ai 2 


1) R. P. GitBert now at the Michigan State University, East Lansing, Michigan. 
2) Department of Mathematics, University of Pittsburgh. 
3) Numbers in brackets refer to References, page 355, 
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Constants 
Lo oil viscosity (centipoise, cp.) 
Mw water viscosity 
P.=P,—P, capillary pressure 
Oo oil density (gm cm) 
0: water density 
k formation permeability (cp. cm® sec gm~*) 
f formation porosity (dimensionless) 


g gravitational constant 
V., velocity at infinity 
n unit normal-vector to the oil-water interface 


t unit tangent-vector to the oil-water interface 
a width of channel 

ee. 
Greg elt 


2. Mathematical Formation of the Problem 


Obtaining an analytic solution to the general problem of hydraulic circula-_ 


tion of two fluids is virtually impossible, since it involves solving the free 


boundary problem in three dimensions. In order to simplify the situation 


somewhat, the following assumptions will be made: 
a) The oil in the oil-cap is immobile (the tilted-water table is stationary). 
b) The flow of water is steady. 


c) The oil-water bearing stratum may be represented approximately as a — 


bent two-dimensional channel. 
d) The fluids are non-compressible. 


e) The generalized form of Darcy’s equation describes the flow of fluids 


through a porous medium. 


Now, regardless of what assumption we have made concerning flow through 
a porous medium, all fluids obey a continuity equation, that is, matter must 
be conserved. We express this by 


Vie V) = s70 %) + $v) + Lev) =—%. (1 


In addition to the continuity equation, we have assumed Darcy’s equation 
must hold. The generalized form of this equation for flow subject to conservative 
‘external’ forces is [4] 


V=—Vd¢, hee d= (6-8), (2) 


a PP nw oe 


Vol, XI, 1960 On a Two-Dimensional Free-Boundary Problem 343 


and ¥ is the scalar potential of the external force field. For the case where the 
external forces are solely gravitational the velocity potential may be written as 


k 
as (Pb+egy), (3) 
assuming, of course, y lies in the vertical direction. When the fluids are incom- 
pressible equation (1) and (2) may be combined to form 


V- (0 V)=—V- (Vg) =-V96=0, (4 


which implies irrotational flow. 
In view of what has been said, we may express the velocity potential for the 
oil and water respectively as 
k 


bo = ial (bo + £004) » (5a) 
and 
bu = <7 bw + 6009) (5) 


If the oil-cap is stationary equation (5a) becomes 


k 
dy = const = —— (fp + £00 4), 
Mo 
or 


bo=—Cory + HoH Po - (6) 


The oil-water interface is characterized by a zero normal flow across it, that is, 
by 
bw _ . 
Oe a (7) 
and a pressure discontinuity, 


Ap le tvestase = Po a Pw fa p : (8) 


By combining Equations (6) and (8) an expression for the water velocity 
potential (at the interface) may be obtained, 


by = : = (04 — Go) ¥ + const , (9) 
WwW 


which is independent of the hydrodynamic water pressure. On the free-bound- 
ary the normal and tangential components of the velocity*) are 


V,=n-V=0, at Vit VG ay te, (10) 


4) Since all future remarks will refer to the water, subscripts will be dropped. 


344 Ratpy J. Bran, Rospert P. GILBERT and Hatcoms KENDALL ZAMP 


where , 
ge AS. 
fap ij (0, 00) , 


VSS 


the x- and y-components may then be computed as 


ax Ge : 
V,=Vie t= Ve ae y cos sin é , | 
and (11) 
da ax \2 : 


Hence it may be readily observed that V, and V, satisfy the following relation- 
ship 
Vee VG eV gp Va a V ees (12) 


We have assumed that the oil-bearing stratum may be represented as a bent 
two-dimensional channel such as in Figure 1. It is worth noting that the 


iy 


Z=X+IY 


Figure 1 


The physical plane. 


orientation of the channel merely introduces a direction parameter in the 
solution to our problem. The angle will in general be obtuse; however, we 
shall first consider the special ‘case y = 2/2 because of the simple formulae 
which result. 

The nature of the complex potential, w= © + 7 VW, corresponding to flow 
in a bent channel, may be seen in Figure 2. The point E’ which is the image of 
E under a conformal mapping of the z-plane onto the W-plane, may be placed 
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on the Y-axis because of the arbitrariness associated with defining a physical 
potential. 


ly 


(Aye NE Se See ade 
be, Ee B Ass 
w-O+iy 
Figure 2 


The complex-potential plane 


To determine the shape of the oil-water interface we shall first obtain its 
differential equation. This may be done by mapping the W-plane onto the 
¢-plane (the plane of the complex potential), and then recalling that they are 


related by 
es dz V.+2V, ’ 
a = = pe ue (13) 


What the image of the physical plane is in the velocity plane depends to a 
certain extent on what boundary conditions are assumed for the free surface at 
B and C. There are certain line-segments, however, whose images are inde- 
pendent of these considerations. For instance, along the straight impermeable 
barrier EDC, the velocity must have a constant direction, hence the image 
E" D"C" (in the ¢-plane) is a straight line in the same direction. If V = co at E 
(which it must be for irrotational flow) the straight line then passes through 


the origin. On CB 
V, 1 
ye y 


(from equation 12), hence C” B” will be parallel to the V,/V? axis; if B isa 
stagnation point B” will be at oo; where C” lies depends on the angle the 
interface makes with the channel wall. The segment £” A" B" is the image of 
the straight impermeable barrier EAB and consequently lies in the horizontal 
direction. It is apparent that E”D"C" and E" A” B" are not in general collinear, 
hence in order for our polygon to close in the ¢-plane B” must be the point at 


_ infinity. This implies B is a stagnation point. 


, 


Ae 


In the following sections we shall investigate the forms our solutions take 
for various boundary conditions at C. As a guide to the solution of the general 
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problem we shall consider first the situation where the channel has a right- 


angle bend. 


3. A Channel with a Right-Angle Bend 
3.1. A Free Boundary without an Inflection Point 


For the case of a channel of constant width, where gy = 2/2, we have the 
following (see Figure 3): 


| 
| 
| eal 


Figure 3 
The complex-velocity plane. 


) At A, V is in the positive x-direction and equals V,,. 
b) At D, V is in the negative y-direction and equals V,,. 
C)LAt Ee, = coat. BV 0, 
d) As we approach B from A, V is in the x-direction, hence BY, must lie 
on the V,/V? axis (since E is at the origin). 


e) On the free-boundary V,/V? = — 1/y, hence the imaginary part of a | 


point on C” BY, is — 1/y (a constant). 


f) As we approach C from D, V is in the y-direction, hence at C V, = 0; 


if V, + 0, then C must be a finite point'). 
This concludes our verification of Figure 3. We now introduce the para- 


meter 
C= cot a = cot1( 5") ; (14) 


in order to describe points on C” BY); it is clear that « is also the angle the tan- — 


gent to the interface makes in the physical plane. 


5) We shall consider the situation Vv, = V,, = 0, a stagnation point at C, in the next section. 


i 
4 


‘ 
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We may use the Schwartz-Christoffel transformation to map that semi- 
infinite strip in the ¢-plane onto the half plane as follows (see Figure 4): 


yi Fad DY” oo Av 
-—_-_----->--——— + OO 
Figure 4 Figure 5 
The t-plane. The T-plane. 
t= —cosh(zwy ¢). 
At pointC”, aylC=—n21t so t= —cosh(— 27) =1; 
At pont E”, xyf=0 so t= -—cosh (0) = —1; 


At point A”, amyf= a so T= —cosh(7,"); 
, Voo pes 


m — BPs ee myt tes uy 
paat J) » wyl= V so T= cosh ( i ) cos| lh 
On E” B”, fis real and positive, so t goes from —1 to — oo. A more convenient 
arrangement of the points A”, D”, E” may be affected by using the cross ratio 
to take 

A"™—>oo, D”—0, and E” +-1. 
One obtains 
t + cos (z y/ Vo) cosh (x y/Vo.) — 1 


ne t+ cosh(xy/Vo)  1— cos (1% y/Voo) (15) 
The points BZ and C” map into 
mM. = cosh (x ¥/ Veo) =i 
a= 1 — cos(% y/Voo) ” 2) 
and 
Co. fs 1+ cos(w y/Vo.) _cosh(x y/Voo) — 1 (17) 


1 — cos(x y/Voo) * “cosh (7,7/Voo) + 1’ 


_ respectively (see Figure 5). The 7-plane (Figure 5) may now be mapped onto 
the W-plane by : 
Voo af 


M4 


W= logT. (18) 


This may be readily verified by considering the points 4%, D'Y and Baa he 
_ purpose of these successive transformations was to map the ¢-plane onto the 
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W-plane, and this may now be accomplished by the transformation 


axa 


= = | 
iS { —cosh(z y ¢) cos(% y/Voo) cosh (7 ¥/Voo) | 
W = log 


|= cosh (x y 6) + cosh (2 y/Voo ). 1 <cosh(ay/Voas 


Our concern in doing this has been to obtain a differential equation for the 
free boundary CB. A point on C” B” in the ¢-plane is described by 


f= + cote, (20) 


1g 


in the t-plane this corresponds to 


t = —cosh(—zi + a cota) = cosh(z cota) , (21) 


and in the W-plane to 


Voo a cosh (x y/V9) — 1  _cosh(a dx/dy) + cos(% y/Voo) \ pany © 
Wome rt ue 1 —cos(ay/Vop) cosh (x dx/dy) + cosh (a ¥/Voo) J’ 22} 
or [since¢zc, = y y + const] 
RES cosh (a dx/dy) + cos(x y/Voo) ) 
a OS hee (x dx/dy) + cosh (x y/Voo) f° (23) 


The constant c can be determined by considering Equation (23) at point B. 


Since dy/dx = 0, cosh(z dx/dy) = co; hence at B, yy+c=0,-and hence © 


c = 0. Equation (23) then becomes 


cosh (x dx/dy) + cos (a y/Voo) 
cosh (a dx/dy) + cosh (% y/Voo) 


= EVV ont (24) 


which is an ordinary, first-order differential equation for the interface. The — 


value of y at the point C can be determined as shown since we have assumed 
there V, = 0, V, + 0, which is equivalent to a zero wetting angle, that is 


dy _ dx ax 
Fister af cosh (x 5] =1; hence 4 
Ved | 1 + cos (x y/V.) oe 

Ye yn 1 + cosh (z y/V,.) 
Equation (24) may now be solved by a single quadrature as 
y 
‘a £ _1[ COS (a-y/Voo) — €~ 772% 00% cosh (2 y] Voo) 

Tonk tate [cosh ‘| rer es | dy. .(26} 


Ve 


3.2. A Free-Boundary with an Inflection Point 


In the last section it was assumed that at C, V,=0 and V, + 0, the alter- 
nate possibility, V,= V, = 0 leads to an interesting situation, Here there will — 


eee ee ee 


ee i 
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be a stagnation point at c, hence the image point C” in the ¢-plane must be at 
infinity. It may then be concluded that the images in the ¢-plane of EDC and 
EAB are the negative, imaginary axis and the positive, real axis respectively 
(see Figure 7). Along the free-boundary BC the condition V/V? = —1/y must 
hold; however By, and CY, are now both at infinity. In order to show that the 
image of BC does not contract to a point it is sufficient to illustrate that the 
image of interior points of the arc BC are not necessarily points at infinity. 
This may be done by again introducing the parameter 


/V../V2 F 
“a= cot 72] == cot (4) : 
y 


it is apparent that «, the slope angle along the interface BC, will in general be 
different from zero unless, of course, no oil-cap exists. It may be concluded then, 
that providing an oil-cap exists, the image of BC folds back on itself at least 
once. Sine there is no physical reason to expect a stagnation point along BC 
other than at B or C there is just one such folding. If F” is the point on the 
image of BC furthest to the left), its inverse image F will correspond to a 
point on BC where « = «,,,,- F is clearly, then an inflection point on the free 
boundary (see Figure 6). 


max* 


De Don 
Figure 6 


The z-plane. 


Now following the same procedure as in the previous section we map the 
¢-plane on the t-plane (see Figure 8) by the Schwartz-Christoffel transforma- 


6) It is assumed that the Re(F”) = 0 in order that the interface BC is single-valued. 
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tion. This may be accomplished by 


Tt 


. 


C=« i (x + 1)-12 c-82 (¢ — t,) dv + B 
er ad T ———— 
=alog|/e? +7 +4 + + + 2a yert+rtB, 
where E” >—1, C4, >0, and B& > co: 


/ v/v? 


V, +/V, 
x 
= 


Figure 7 Figure 8 
The complex-velocity plane. The t-plane. 


A listing of how certain points map from the t-plane to the ¢-plane may 
be found below: 


Points in t-plane ¢ as computed From ¢-plane 

BY: (t = +00) t=alogoo+2at, +P By, = 00 =i . 

By: (t = —oo) € =alog(—oo) + 2at, + 8 Bi = 0° 
=alogcootani+2at,4+ 8 ' 

C": (c=048) C=alog> +2at,+8 Ch, = 00 i 

1 ‘ 
C™: (t= 0— 2) C=alog5 —2ait,co +f Cl =—ocoi 
mM 1 
E": (c= —1) ¢ = alog(—3) +6 E"=0 


th ; 
=alog>+ani+ B 
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Some of the arbitrary parameters may be computed now by comparing various 
terms. At E", ¢=0, hence 6 = ~a« log1/2 — «a7. By comparing BY, f = oo, 
with € computed from t = — oo, one has 


€ =co=co+ (ami) + (20 t) —alog> — asi; 


hence t, is real. Comparing B”, € = co —i/y, with ¢ as computed from 
t = + oo, one has 


ee 


= i : 
¢ =0O —— 1=0ce+2at, —a log, —ant; 


, 


consequently 


me a a 1 i 
Carre and fee oe Oe eae 


In the same manner, at C4, 


Peay ie eco e siege 
Z eee oot ih a oe ne Da 
hence t; > 0. Again, at C4, 
4 
= 
= OO —-— 
c oe 
i 1 1 2 1 i 
fore) = : ae +2 tT, co + —— log2 — —; 
yp ay ©2 a ae ny © y 


hence, as concluded above, t; > 0 (a check). The point’) ¢; may be calculated 
from t; as follows 


= Ho. 1 Ti oS, 1 Tt ——— 
l= — a+ as [log +etr+ sl} 4+ 2 yet x + log2]. (26) 

As before, the t-plane is mapped onto the T-plane, which may then be mapped 
onto the plane of the complex potential by 


Von & 
Vem log. 


_ First use the cross ratio 
GE — op 14+, 


ere feoe i 1+ Tp’ 


_ which takes tp > 0, co > 1, and ty > o; it follows that 


7) We shall adopt the notation of using letters as subscripts to indicate image points in the 
- remainder of this section. 


‘ir ty ly ct Dele tok Ne 
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Now, on the interface CFB, 


Me py o> aie tie 2 
BE oe ne ae y I y dy’ 

and 

W=@O+ yy+ constant; 
consequently 

~ Ves @ t—TtD) K 
I y IU log 4H + TA 7 mee 

At 


Ay Mi Vinge 0 sivas ee co)= 0; 


therefore K = 0, and we have 


y 


a ty) tA elt IV 994) 


1 — elt7/Vooa)y 


ZAMP 


(27) 


Using equations (26) and (27) we may eliminate t, and thereby obtain a dif- 
ferential equation for the free-boundary; this equation can not be solved 
explicitly, however, because of the undetermined parameter t;. Consequently, 
we must show whether C is a stagnation point or not, and if C is a stagnation — 
point, t; must be determined. From physical arguments it may be realized, 
that if the free-boundary has an inflection point at F this corresponds to an 
unstable equilibrium; we conclude that a stagnation point can not occur at C. ~ 


4. The General Problem 


In the preceeding sections we have considered y = 2/2; that was done — 
primarily because of the simpler formulae which arise for that case. Now, we 
shall assume, as in section 3.1, that C is not a stagnation point (V, = 0, V, + 0); 
consequently, we obtain the complex velocity plane as shown in Figure 10. 


Figure 10 may be verified by considering the following: 


a) Along EDC (a straight impermeable barrier) V has constant direction, 
at E; V = oo, consequently the image in the ¢-plane is a straight line through 


the origin. 


b) On CB (the free boundary) V,/V? = —1/y, V, + 0, hence C" is a finite 


point [1/y cot (x — gy) —1/y], and B" is the point 


(S52, ea 9), 


The shear-transformation 


E=€ — cot (x — 9) Im (¢) = Ve +1 Vy a cot (a — g) Vy 


Ve V2 


takes By, E" D’ C" BY, (Figure 10) into a semi-infinite rectangle (Figure 3). 
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Ww 
Ga 
WwW 


Figure 9 
The T-plane. 


i Wy 
y2 


Figure 10 


The complex-velocity plane. 


We may now proceed as in Section 3.1; 
t = —cosh(zy é) 


maps the semi-infinite rectangle on the upper half of the t-plane, where 


and 


The cross ratio is again used for a more convenient locating of Tg, Tp, and T48),; 


T= (= £OS[(% 7/Voo) sin (x — ¢)]) [cosh (7 ¥/Voo) — 1] 
(t+ cosha y/V,.) (1 — cos[(a y/Voo) sin (7 — ¢)]) 


(See Figure 5) where 
. T= cosh (zt y/Vo) — 1 
B 


1 — cos [(a y/Voo) sin (a — @)] 


1 + cos[(% y/Voo) sin(% + @)] | cosh (z y/V..) — 1 
1 — cos[(x y/Voo) sin(a — y)]  cosh(z y/V,.) + 1° 


To = 


o 


8) If p = m/2, this equation reduces to the one in Section 3.1. 
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On the interface ¢ may be given as 
it il 
fF = —— 4+ — cota, 
° ¥ rs Me 
where « is the parameter previously introduced and « varies from a — @ to 0; 


consequently, 


t = —cosh(xy é) = —cosh[w y € — wy (Im¢) cot (x — 9)] 


Poteet V, il _ ae 
(Im< = yr eos, and cota = a ; 
t = —cosh[— zi +acota+a cot (a — 9)] 


= cosh[z cota + 2 cot (x — g)], 


= cosh (z cot«) cosh[z cot ( — ¢)] + sinh (z cota) sinh[z cot (a — g)] . 
The combined effect of the transformation on the interface is 


Veo & 
Up ee prege Clee! 


Vana 1 { (4, cosh (x dx/dy) + K, sinh (a dx/dy) + Ky] (K,— 1) 


Die ay 08 \ 0c crab (x dx/dy) + K, sinh(x dx/dy) + K,] (1 — it 7 oe 


K, = cosh[{z cot (a — q)], K,=sinh{[z cot (x — g)], 


Ks = cos [eo sin (7 — )| ou Sandy cosh (Fo . 
Recalling that ¢z- = Ay + constant, and using the condition dy/dx = 0 at B, 
one may obtain the following relationship: 


Voo @ { K, cosh (a dx/dy) + K, sinh (a dx/dy) + aa 


Vil Os NOB ae acta dx/dy) + K, sinh (x dx/dy) + Ky | 


Using the identity 
K, cosh (= a + Ky sinh (= =\= ae re {cosh (x= ae + tanh? ie }t. | 


where Kj = K3, one may solve equation (30) by a single quadrature as 


y 
1 as RYI[V oo 4 
eae t= | |cosh3/ : oe ele By dy 
A A Ke — i k3 et VIVog4 aed 
Ve 


i (tanh— Zz) (Y = yee {am 


Ee gs 


_ 


‘4 Y 
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The boundary condition y=y, at x= *, may be cetermined from equa- 
tion (31) by setting dx/dy = cot (p — a) at point C. One obtains 


Ven Ke — Ke Ke 
ee hg ee ~ ey and x,=y,cot(p—z). (32) 


Conclusion 


The problem of determining the free boundary between two fluids, in a 
bent channel, under dynamic conditions has been investigated. For a channel 
with a right-angle bend the possibility of the interface having an inflection 
point was also considered. This situation appeared to correspond to an unstable 
equilibrium, hence when the general problem (a channel with an arbitrary-angle 
bend) was attempted it was assumed that an inflection point could not occur. 
An ordinary, first-order, differential equation for the free-boundary was ob- 
tained, which was pointed out could be solved by a single quadrature. 
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Zusammenfassung 


Untersucht wird, wie in einem Kanal mit Abbiegung unter dynamischen 
Bedingungen die freien Randwerte zwischen zwei Flissigkeiten AL bestimmen 
seien. Fiir einen Kanal mit rechtwinkliger Abbiegung wird insbesondere die Mog- 
lichkeit in Erwagung gezogen, dass die Zwischenflache einen Wendepunkt besitzt. 
Dieser Sachverhalt bezieht sich auf einen unstabilen Gleichgewichtszustand ; fiir 
den allgemeinen Fall (Kanale mit Abbiegungen in beliebigem Winkel) wird deshalb 
angenommen, Wendepunkte kénnten nicht vorkommen. Gewonnen wird eine 
gewohnliche Differentialgleichung erster Ordnung, fiir die sich herausstellt, dass 
sie durch eine einfache Quadratur gelést werden kann. 
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A Type of Solution for Thin Shells in the Form 
of an Hyperbolic Paraboloid 


By Girpert H. Becuin, Genéve 


Introduction 


Today some interest is being shown in the shallow shell theory (2, 4, Sifu 
One interesting problem is that connected with the hyperbolic paraboloid, for 
this surface has a simple algebraic expression, and it is a surface of negative 
total curvature. 

A system of differential equations which describe the state of stress and 
deformation of shallow shells has been indicated by MARGUERRE [3]; this 
system takes membrane, bending and non-linear effects into account. Russian 
authors [1, 4] have derived similar systems of differential equations for the 
linear case. ; 

We recall that by ‘shallow shell’ is meant a shell whose middle surface ~ 
2(x, y) is such that the quantities 


cS Sey (ae 


are negligible compared to unity over the entire region occupied by the shell. 
In practical applications the use of this ‘shallow shell’ theory is allowable as 


long as 
(a5): 


depending upon the purpose of the analysis, and its required accuracy. 


Statement of the Problem 


For convenience we shall use the notation and the equations of shallow — 
shell theory as presented by E. REISsSNER in his article SHE 


Oxy is a Cartesian system of coordinates, its origin 0 lies on the middle 
surface S of the shell; S is described by the equation 


a= 22, y= ay. (1) 


i) Numbers in brackets refer to References, page 367. 
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This is a particular type of hyperbolic paraboloid, for its two directing planes 
are perpendicular. The z-axis is parallel to the intersection of the directing 
planes, and the generators of the surface (1) are projected onto the x y-plane 
as the lines x = constant, and y = constant. The lines of curvature project 
onto the same plane as the curves: 


u + v=constant, with sinhu = °~ sinhy = ©» 
ab? wi lo 

AZ 

-b | 
Ay 
Figure | 
Middle surface of the shell. 
We consider the portion of the surface (1) bounded by 
Ca Hwee, and v0, pe =o. (2) 


The edges of this region are generators of the surface, and straight lines 
(Figure 1). This surface we take to be the middle surface of a thin shell with 


uniform thickness h. 
The conditions of equilibrium of a small element yield 5 equations for the 


force and couple resultants: 


ey fgg 2 Nn) tangential stress resultants, 
OF, O53 transverse stress resultants, 
MM ys Me, = My: flexural and torsional couples. 


‘The load is defined by its components in the x, y, and z-directions: 9,, q,, -- 

The force-displacement relationships in classical theory (Love’s first ap- 
_ proximation) yield 6 additional equations. All these equations can be reduced 
to a system of 2 simultaneous differential equations for the Airy stress function 
F and for the displacement component W. 
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In the case of the linear theory these 2 basic differential equations read as 
follows (q, = 4, = 9): 


2G 2 
D-P2V2Ww oS a V2 V2 Pe pe aaa Ee 


= —~-. . (3 
= 42~ Gb Ox Oy’ ab 0x oy (3) 


The problem is to find a set of 2 furctions F(x, y) and W(x, y) satisfying (3) and © 
some boundary conditions to be prescribed. 
Remarks 
It would be convenient to get rid of the term q,. Thus we set: 
W=W,+60, F=R Yo: (4) 
and substitute in (3). We note the following: 


a) Whenever a function /(x, y) can be found which satisfies both the 
equation of linear membrane theory (g,: given loading): 


BOO 0 (5) 


and the condition: 
272 =0; (6) 


then the system (3) can be reduced to the following: » 


oars :* 26 ODT e O2@ 
Uns fina 1 Ox Oy ’ 
PPeO=Eh dale ead "4 
= “Ox Oy 2 “Ox Oy * 4 
The parameters k, and ky are given by the following expressions: { 
; 
4S ke 74 (0 : E hs 4 
ky ey Ra = an Eh, with D= Ge (Ry, Ra, D> 0), (8) 
b) Whenever a function W,(x, y) can ‘be found which satisfies both the 
equation of classical plate theory: 
D-V?02W,=4,, (0) 
and the condition | 
OW, 
Ox Oy 04 


then the system (3) can be reduced to the following: 


— 
— 
c=) 
— 
RS LW ee ee 


02F 
Me Ea le a ae NE eM Sets ie 


— 
— 
Ss 
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From now on we assume the load to be such that one of the above reductions 
is possible. 

Within the assumptions made above the problem of finding the state of 
stress and deformation of an hyperbolic paraboloid subjected to the load q, 
— and supported in a way still to be prescribed — amounts to finding a set of 


_ 2 functions, say W(x, y) and @(x, y), or 0(x, y) and F(x, y), which satisfy the 


system (7) or (11), and appropriate boundary conditions. In the latter the 
contribution of /, or 1V,, must not be forgotten. 


Solutions (general considerations) 


Let us assume that a solution of the system (7) exists. For the time being 
we disregard the boundary conditions. We call this solution T(x, y): 


T = [W(x, y); P(x, y)], 


T is a vector whose components are W and @. 


Property 1: if ie and Ty are solutions of the system (7), then fi + A in 
is also a solution of the same system (A = constant). 


Property 2: if W and @ are the components of a vector T, which is solution 
of (7): T= (W, @], then 


(W+k,?@d) and ®—k, 2 W) 
are also the components of a vector Ve 
T, =((W+k, 28); G—k, 2? W)], 
which is also a solution of (7). The arbitrary parameter / has the dimension of 


a length. 
The operator //? reduces to the following: 


72 = ed ioe (12) 


because of the assumption of shallowness and of the choice of coordinates. 
Simple solutions in terms of polynomials are obtained by direct substitution 


_ and identification of coefficients. One type of solution reads as follows: 


_ 


- 


~~, 


@ = y [a, %° + P,(x)] — x(a, ye Q,(9)] + Pa(x) + Qal(y) . 
W = y [c, x5 + U,(x)] — xfer v5 + Valy)] + Ug(x) + Voly) - 


P;(x), Q;(v), U;(x), Vj(y) are polynomials of degree 7 in % and in y. The poly- 


nomials in the above expressions are not entirely arbitrary, some of their 


- coefficients are interrelated. This type of solution is of limited interest. 
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A Solution of the Levy Type 


In the plate problem, in classical theory it is customary for the case of a 
rectangular plate with 2 opposite edges simply supported to take, as shown by 
M. Levy, the solution in the form: 


[o,e) 


W= Dy. Y(n, y) sin ane : 


n=1 


For the problem of the hyperbolic paraboloid under consideration, in linear 
shell theory we take a solution in the form: 


W = w(y, a) sina x, 
= - » (13) 
| © = x fly, @) cos «x; with #27 =F h D’. | 


We substitute these expressions in (7) and obtain 2 simultaneous ordinary 
differential equations for the unknown functions f(y) and w(y). These equations 
read as follows: 


df a2} dw . 
i dtw 2 Vw ~ af 
Lw) = Gar — 20 Ge totw= haz, (15) 
where the constant k is a positive number: 


The operator L( ) applies both to f(y) and w(y). In order to obtain a solution 
for this system (14, 15) one can eliminate one of the unknown functions, say 
w(y), and obtain an eighth order ordinary differential equation for the remaining 


function, say /(y). However, we prefer to proceed as follows: 
let us put 


fy) = a; ens wy) = b; ee (17) 


and substitute in the system (14, 15). This yields 2 simultaneous relationships 


which read: 


a (M — 2a A tyot\ eng (72 ayr= db kad, (18) 
b; (AA — 20? 12 + a4) = b, (2? — 02)? =a kad. (19) 
There are 2 alternative possibilities: 


a,=b,, and a,=—b,. 


‘ 


——— 


OO — 


PN 
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If we choose a; = 6,, then the corresponding characteristic equation reads: 
(A? — o?)2—kadA=0, casel. 
If we choose a; = — b;, then the corresponding characteristic equation reads: 
(A? — a)? + kad=O0, case 2. 
We introduce 2 non-dimensional parameters € and f: 
A=&a, k=fBo?, freal>0: 


and write the characteristic equations as follows: 


(A=—3)? = Be= 0. (20) 


~ 


In case 1 there are 2 real roots &,, €,, and 2 complex conjugate roots (—&,), 
(— &,), with negative real parts. 

In case 2 there are 2 real roots (— &,), (— &), and 2 complex conjugate 
roots &, €,, with positive real parts. 

It is possible to establish the following inequalities: 


ae, 1+e, and §=—1-— pz, 
thee Opel; 0O< pee tp) <enu; 
Hd f=yte0, and &=y—206, 
then 0<8; 0<#< 56; 1<y. 
In case 1 the following expression: 


f(y) = a, e®” + aye” + a, e-* cosa y + a,e-** sing y 


(21) 
w,(y)=fh(y), where t=ay, o=a0, 
is a solution of the system (14, 15). 
In case 2 the solution reads as follows: 
- haly )=a,e-*” + aye” + age” cosay + age" sing y, we(v) = —fe(y). (22) 


meh the postulate 


W =wl(y,a)sinax, B=xfly,a) cosax, 
‘the general solution of the system (14, 15) is written 


Hy) =hly) + faly), wy) = wily) + woly) - 
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The above expressions can be written in a fashion that brings out their odd 


and even parts: 


H(y) = a, sinhA, y + ag sinhd, y + a; coso y sinht y + 4, sino y cosht y 


+ ay coshA, y + a4 coshd, y + ag cosa y cosht y + ag sing y sinht y; 


w(y) = a, sinhd, y + a, sinhd, y — ag cosa y sinht y — ag sina y cosh tT y 


+ a,coshd, y + a3 coshdy y — a, cosa y cosht y — a, sino ysinhty . 


The coefficients a; are determined by the boundary conditions. 


Another Solution of the same Type 
The postulate (13) may be replaced by the following: 


W = w(y, «) cosa x, 
(i (24) 
|o= Hy) oa) sina x? with x*= & WD ; 


This leads to a system of differential equations for the unknown functions — 


w(y) and f(y), which reads: 


Ge a 
Lif) = 4 =e ems + at f = ha, 

d*w a af a 
L(w) = Fa -2 2 dy? te 


f 
U 


In a manner similar to the one used above, it is possible to derive for f(y) and ~ 
.: 


w(y) the following expressions: 


Hy) = by eh” + by e“*® + By c®” 4 bye“? 


+ 6; e"” cosa y + bee’ sing y + b, 6-7 cose 4) +- 6, e “sing y; 


— 
dD 
(=) 

er 
a 


w(y) = —b, e*” + Bye“? — B, eh? 4 Oye 


ty pees 
+ bs; e*” cosa'y + bg e*’ sing y — b,e-** cosoy — bge*’ sina y. 


ie so aM 4; =a&,t=ay,o=«6 are exactly the same aes a 
ah he system (25) leads to a characteristic equation identical to (20); i 
eads 


(6? = 1) 82h Bete OF 


a _ 
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For purposes of application we write the functions f(y) and w(y) in a way that 
shows their odd and even parts: 


f(y) = 6, sinhd, y + 6, sinhA, y + b; cosa y sinht y + b, sino v cosht y 
+ 6, coshd; y + by cosh, y + bg cosa y cosht y + bg sino y sinht y ; | 

w(y) = — by sinhd, y — by sinhA, y + bg cosa y sinhr y + bg sina y cosht y oe 
— b, cosh, y — 6; coshd, y + b; cosa y cosht y + b, sino y sinht y. | 


The coefficients b; are determined by the boundary conditions. 


An Immediate Extension 


In our postulates (13) and (24), the parameter « was left arbitrary. A genera- 
lisation of (13) presents itself immediately in which the solution can be written 


as follows: 


W = dD, wy; n) sin=2* ; 
n=1 


(rd 
| B= Dd’ x fly, n) cos ===, 


n=1 


(28) 


ns 


The same holds true for (24). 
It must be clearly understood that for each value of the index m there are 
4 roots &,, &, &3, &4, of the characteristic equation (20), and in general 8 terms 
in the expressions for f(y) and w(y). The tabulation of the roots €; is an involved 
task, for their magnitude depends on the geometry (form and thickness) of 
- the shell under consideration. More precisely the magnitude of the roots &; is 
governed by the parameter f which reads: 


eee br (455) i ; (29) 


B ies: 7 


We note also that the boundary values of the various quantities (couples, 
; displacements, etc.) along the edges y = constant appear as ordinary trigono- 

metric series. They can be handled in the usual fashion. However, along the 
edges x = constant, the boundary values as derived from (28) involve series 
of terms of the type (23). It is not known if an arbitrary function can be 
_ represented over (— 8, b) with the help of such terms. In specific cases one may 
_ have to seek an approximate expansion. 


~ a» te 
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Application 


Consider the portion of the hyperbolic paraboloid 


ah Sous 


limited by the planes 


(Figure 2) 


a 
*=0, =a, and y= 2? 


“vy 


1 


Figure 2 
Portion of hyperbolic paraboloid. 


With this choice the parameter f may be written 


piel fon cea) 


ZAMP 


(30) 


For the sake of simplicity we take the factor in the brackets to be an integer, 


say 36. This corresponds to a shell thickness of about 2 = 0-009 65 a. 


Then the roots of the following equation 


36 
(Ph) ee 


must be determined. 
Computations show these roots to have the following values: 


= lie 


&,= 1+ ¢ = 3-495 357...; &, = 1 a= 0027735 0res 
&=y +26 = 1-761 546... + 7 2°685 556...; E,=y—10. 
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is 2: 
€, = 2-370 409-..; & = 0-108 509 9...; 


g=y +76 = 1-439 459... + 71-347510...; &=y—i6. 


&; = 1-946 965...; €,= 0-225 2704...; 
&s = y +16 = 1-086 117... + 7 1-048 981...; &, =y —id. 
The shell is assumed uniformly loaded: 
Ge = WY =, 9, = p = constant. 


Along the edges x = 0 and x = a, we choose to impose the following boundary 
conditions : 


“=0 and =a: W=0, M,,=0, i.e.simple support: | 


(31) 
u=0, N,,=0; | 
i.e. edge stiff in the x-direction. 
The appropriate form of the solution is as follows: 
| W = D*wly, «) sina x + W(x, 9); | 
= (32) 
|@ = x Sj (y, a) cosa x, with: 2° == 2D. | 


The conditions (31) are automatically satisfied provided we take 


a=", (0,9) = Mla, 9) = 0, 


and atx = Oand*x=a: 
ow, | at = 10 
eee PI al a ee es 


The function W,(x, y) must also satisfy the conditions (9) and (10); it is then 
a function of x alone which reads: 


4pat*wl NT x 


W(x, 9) = W(x) = =P @t — 248 + ax) = De sin, (33) 


24D 


where is an odd integer: = 1, 3, 5,.... 
For w(y,«) we take an even function of y in order to obtain the same 


- boundary values on both edges y = + 4/2: 
w(y, «) = a, coshd, y + a, coshAy y + + +5 (34) 
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the corresponding general expression for /(y, a) reads: 
f(y, «) = a, sinhdA, y + agsinhA,y ++". (34') 


The arbitrary constants a, will be determined with the help of the boundary 
conditions at y = + a/2, still to be prescribed. These conditions involve stresses, 
couples, and displacements. For these quantities we have the following general 


expressions : 
N,, = * D>, (—22) fly, a) cosax, Nyy = x Sa f(y, a) sinw x, 
M,,=—D “2 ae a) — va? w(y, x) sina « + » W,'(x)], 
M,,= —(1—») D Sia w'(y, «) COSa x, 


4 


b= a ox = ; (i, a) + va? f(y, a) sing * — —- ZY W(x, y) + const , 


x 


Ta yy (-#f 10 aldy—yvri ly, x) cosa % — Se x W(x, y) +1(x) , 


where uniform convergence of the trigonometric series has been assumed. 
The arbitrary function /(x) is related to W,(x) in the following manner: 


I(x) = 26. [me dx. (35) 


We are now in a position to impose at y = + a/2 boundary conditions compat- 


ible with the postulates (32) and (34). We may require e.g. these edges to be 
simply supported: 


W (x, £5) = \' (+5, x) sina x + W,(x) =0, 


Jide (x, +5) =0, or > 2" (7, «) S053 = (0) - 5 


We may also require these edges to be stiff in the y-direction: 


v (x, +5) =0 and Ns (x, +5, = 


These 2 requirements lead to the following equations: 


a 


al2 


ps lect oe ) dy — vf (+ - Beem eo 


Sat (£5, a) sina =0, 


—_ sa. 
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As an approximation we take only the first harmonic (n = 1). Thus we get 
4 equations in 4 unknowns for the coefficients entering the solution (32). Com- 
putations show these coefficients to have the following values: 


4pa’ , 
a; = pe a; , with 
ay Sena hOs5 3G ay = —0:997 519...: 


a, = —10-8 -0-2778...; a, = 10-3- 0-714 5... 


The expressions for the stress resultants and the couples are then easily 
written. 

Remark: It is seen that the boundary-conditions as imposed in the above 
example lead to a state of stress which is unlike that of a flat plate. 
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Résumé 


Pour le probléme du paraboloide hyperbolique une solution est présentée qui 
se base sur la théorie linéaire des voiles minces peu incurvés. La solution prend en 
considération l’effet de la flexion; c’est une extension du type de solution indiqué 
par M. Lévy pour la plaque rectangulaire (théorie classique). Un exemple d’appli- 
cation termine cette étude. 
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The Minimum Weight Design of Circular Sandwich Plates 


By Greorce Eason, Newcastle upon Tyne, Great Britain*) 


1. Introduction 


The minimum weight design of circular plates was first discussed by Hop- 
KINS and PRAGER [1]?) for the Tresca yield condition. FREIBERGER and TEKIN- 
ALP [2] assumed the von Mises yield condition to obtain the solution for a 
simply supported, fully loaded plate. Design for minimum weight was first 
made systematic by DRUCKER and SHIELD [3, 4]. The direct design procedures 
developed by them were used in [5] to solve the problem of the minimum 
weight design of sandwich and solid circular plates assuming the Tresca yield 
condition and the associated flow rule. 

The direct methods due to DRUCKER and SHIELD are used here to investigate 
the minimum weight design of a circular sandwich plate when the von Mises 
yield condition and its associated flow rule are assumed. The material of the 


plate is assumed to be perfectly plastic, non-hardening and uniform so that the © 


design for minimum weight is the same as that for minimum volume. The 
minimum weight design is to be such that the whole plate is about to collapse. 
The plate is assumed to be simply supported or built-in at its edge and is 
uniformly loaded over an axially symmetric area. The results obtained for the 
case of simple support are identical with those for the Tresca yield condition. 


When the plate is built-in the results are more complicated and appear as © 


integrals which have to be evaluated numerically. Comparison is made with the 
results for the Tresca yield condition when the plate is fully loaded and built-in. 
The results for the two yield conditions are very similar in general nature and 
in view of this the much simpler results of the Tresca yield condition are likely 
to be adequate in a practical design. 


2. The Basic Equations 


The plate is assumed to be horizontal and to have radius R. With radial 


co-ordinate 7 the equations of equilibrium when the applied load has axial 
symmetry are 


“7 Q)+rp=0, (1) 
- a VM)-N=rQ, (2) 


1) King’s College, Department of Mathematics. 
*) Numbers in brackets refer to References, page 375. 


eo 
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where M is the radial bending moment, N the circumferential bending moment, 


Q the shear force and # the applied luad, considered positive when it acts 
downwards. 


The von Mises yield condition for plates is 


f=M?—MN+N®2- Mz=0, (3) 


where M, is the fully plastic bending moment. Associated with the bending 
moments M and N are the curvature rates x and 4 given by 


dw of 
sera ogee WON) 4) 
dw of 


Sade ta NSM), (5) 
where yw is a non-negative parameter and w is the downward velocity of the plate. 
Equations (4) and (5) may be combined to give 

aw 


oe gnu - 1.28 eum) 6 


The sandwich plate considered here has a central core of constant thickness 
H and two identical face sheets of thickness ¢ («<< H). The central core carries 
no bending stress so that the fully plastic bending moment is given by 


M=o0,Ht, (7) 


_ where oy is the yield stress of the face sheets. We wish to determine the mini- 


mum volume of the face sheets, allowing the thickness ¢ to vary. The rate of 


: dissipation of energy per unit volume of the face sheets of the plate due to 
_ plastic action is given by 


=~) | COes, 


oh) ie ta 


Ra Se RRR 


18) Mzu+t+NaA 
t = (8) 


where D is the rate of dissipation of energy per unit area of the middle surface. 
It has been shown in [2, 3,4] that D/t is constant for a minimum volume 


design. 
The yield condition (3) is satisfied by the relations 
eS eee ge IME Ue es 
M = — = sing — 8), ciety sin( «+ 6), (9) 


where the parameter 9 has the value 2/2 when M=N=M,. Using the 
relations (9) and (7), equation (2) now becomes 
dé 


vt a cos ( 


1 CL ed a an ai _ r@y3 
6% 0) —7 7 sin( 5 x 0) t cos 0 a (10) 
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and equation (6) takes the form 


; ie 1 
a cos( 7 — 6) ay <* cos(_ a + 6) == (() (11) 
Equation (8) results in 
2 1 he rm A 3 
a sin( 5 It 0) ; — sin (-¢ m4 6) Je (12) 


where A is an undetermined constant. Equations (11) and (12) may be combined 
to give 


dw 1 

7 = A cos (5 z+ 6) 5 (13) 
1 dw 1 
; eee cos (_¢ % — 6), (14) 


and eliminating w between these equations results in 


y ae sin (= mu — 0) + /3 cos 0 =. (15) 


3. Solution of the Basic Equations 


In this section solutions of equations (10), (14) and (15) will be obtained. — 
The procedure is to solve equation (15) for 7 in terms of 0 so that equations (10) 
and (14) may then be solved for ¢ and w respectively. Two solutions of equation | 
(15) are required here, and these will now be considered. 

(i) Particular solution. 

One solution of equation (15) is obtained by assuming that 6 is constant 
so that 


Gs ee (16) 
Equation (14) now gives 
1 
ecaboateode qe * (17)9 


where B is a constant. An expression for ¢ is obtained from equation (10) in the 
form 


om Ht=/Qdr+C, (18) 


where C is a constant and Q is determined by equation (1). 
(ii) General solution. 


In general of course @ is not constant and writing equation (15) in the form 


2 ar a 
Pe do Rk is 


regent neem 
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results in 
= Dsee6 exp(— a fe (19) 
beeen 
where D is a constant. Using equation (15) it is possible to rewrite equation (14) 
in the form 
dw _ A Dcos(a/6 — 6) sin (77/6 — 8) 2 alV3 


dO V3 cos?0 
so that 


w=E-—AD (2 /3e"V5 4 tan 6 e-3 + 1/3 [tang e-AV3 a6}, (20) 


where F is a constant. Similarly equation (10) may be written in the form 


dt 1 1 vQ 
dé - ts 5 cot (5 a 0)| ~ 26H cosé ’ 
so that 
oo H tsin(— % — 6) ee ae : [7 Qsin( 5 — 0) sec 6 e V3 gg , (21) 


“ 


where F is a constant. 


aa Particular Problems 


‘The general results obtained in the previous section will now be applied to 


' some particular cases of interest. 


(i) Simply supported plate, central loading. 
The problem to be considered is the minimum weight design of a simply 


supported circular plate of radius R which is loaded uniformly over a central 
~ region of radius a with the remainder of the plate unloaded. Thus # is constant 


for 0 <r< aand is zero for a< 7 < R. The conditions to be satisfied at 
the edge of the plate are M = 0 and w = 0. The quantities Q, M, w and dw/dr 
must be continuous at r = a. The condition M = 0 at r= R is satisfied by 


taking ¢ = 0 at the edge of the plate. The solution of this problem is obtained 
by assuming 6 = 2/2 everywhere. Applying the conditions given above to 


oo? ee, ee: 


Ay 


Pe eee IE ACA 


equations (1), (17) and (18) leads in a straightforward manner to the results 


6=in, 0S7SR, 
w= A(R—A), OSTSR, 
1 
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1 1 R 
Pe 1) a Pe hoe NSPS @, 


o,Ht=; 
| 3 pattog®, ay ie 


These expressions are identical with those obtained by Onar et al. [5] 
assuming the Tresca yield condition and the associated flow rule. The results 
for the special case of the fully loaded plate (a = R) were previously obtained 
by FREIBERGER and TEKINALP [2]. 


(ii) Built-in plate, central loading. 

We now consider the problem of a built-in circular plate of radius R which 
is loaded uniformly over a central region of radius a with the remainder of the 
plate unloaded. The conditions at the edge of the plate are w = 0 and dw/dr = 0. 
The second of these conditions implies that the edge of the plate corresponds 
to the point on the yield surface at which M has its minimum value and 
6 = — a3. 

The solution of the problem is obtained by assuming that 6= 2/2 in a 
central region of the plate of radius } and that outside this region 6 = — w/3. 
The bending moment M must be continuous at 7 = b, and this condition is — 
satisfied by taking ¢ = 0 at y = 0} and allowing 6 to be discontinuous tkere. It 
was shown in [5] that dw/dy must be continuous at r = b. These conditions 
are all satisfied by allowing 0 to have the value 7/2 as 6 is approached from 
below and to have the value — 7/6 as b is approached from above, so that 
from equation (14) 


The quantities Q and w must be continuous at r= } and also M, Q, w and 
dw/dy are continuous at 7 = a. 
Using the condition y = R when 6 = — 7/3 equation (19) gives 


= Ri sec exp(——a), Direee aoe iy (22) 
where 
$= 5 ex Z 
2 P/ 3 | , ' 
Since y = b when 6 = —z/6 it is found that 
B75 Pash | 


so that 6/ = 0-653 in contrast with the value of 2/3 obtained in [5] for the 
T'resca yield condition. 


| 
| 
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The remaining analysis is straightforward and it is found that equations (1) 
and (17) in the region 0 <y < #4, and equations (1) and (20) in the region 
b <r <R lead to the results 
ee io) 4 
4 4/3 
if 


| af3V3 _ a6V3 __ ay it \ Ms 
Be Se 31, (—% Mp O<r<b 
ae 


mt 
ie As R? 3 ers V3 _ 2y3 e~IV3 _ tan 6 e~ V8 _ 3 1,(0)\, D7 SRS 


1 , 
exes (Oe ee 
Oa 
1 2 
|< pa, OP SIRE 


for the deflection and shear force, after writing 


For a < 5, equations (18) and (21) give 


4 p(a@—r) + 5 patlog?, 0=75¢@, 


¢ a b 
; — p a®log— iy ee 
o,f t= asi ae asrse, 


. a 


par V3 eMll® [ —olVe x glo Ve ee 1,(0)! 
Fsin(aJ6 — 6) 1" é 1, ( 6 zr) ae 1( My 


mV 


For a = d equations (18) and (21) lead to 


{bb (R—), 0SrSd, 


p R25 ellV3 we ghee 
Facer Ne PSS 


paryael® | os _ ,-alVa_ 7 1,(6 
sim (a/o= 0) le F ia) + 16) 


eke PA a Et 


p Res cll8 <rv<R 
+o mae ee 


where J,(0) is defined by ; 


© sin(n/6 — 6) e295 ag 
cos? 6 


3 
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and § = a when r = a. Equation (22) gives a relationship between a and « of 
the form : 
a* = 5 R* seca exp(- a i 
The variation of o, H t/p R? with 7/R for both the Tresca and von Mises 
yield conditions is shown in the table for the case in which the whole plate is 
loaded. 
The volume of material in the face sheets is given by 


R 
Va4afrtdr, 


0 


so that for a fully loaded built-in plate 


Se D Q 

pre! 2ap Rts? sin (1/6 — 6) cos(x/6 — 9) e248 gg 
: V3 cos? 0 

— 7/3 


=O-lll ap RS, 


in contrast with the value of 0-117 2 p R* for o) H V obtained in [5] for the - 
Tresca yield condition. The case of a point force of magnitude P is obtained — 
as a limiting case of the above by letting 2 pa? > P as a >0. The volume 
for this case is found to be 


i 1 = ‘3 5 > iL 
o, HV = 3 Pb+ . P R*s/3 ee = ero V3 _ I, (-;>)| 
= 0:324 PR? 


in contrast with the value of PR?/3 obtained in [5] for the Tresca yield condi- 
tion. 

It remains to be shown that the parameter “ appearing in equations (4) and | 
(5) is non-negative. In the problems considered here the parameter 6 lies 
between — 2/3 and x/2 so that cos (z/6 — 0) = 0. Hence, the penis in 
appearing in equation (5) will be non-negative if the constant A appearing 
in equation (14) is non-negative. In all the problems considered here, A has_ 
the same sign as the velocity, which due to the conditions of loading will be 


positive. Since w is non-negative the velocity fields are admissible velocity 
fields. 
(iii) Annular loading. | 
Results for a plate which is simply supported or built-in, and loaded over 
the region a <r < R only, may be obtained by methods similar to those used 
above. As for the Tresca yield condition (see [5}), the thickness ¢ is constant in 


24 
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the region 0 Sr <a for both methods of support and in particular if a > b 
the thickness is zero in the region 0 <= x < a for built-in support. ve 


The variation of o, H t/p R® with r/R for uniform loading (a = R). 


0 Ht/p R? O, Ht] p R? 
r/R for the Tresca for the von Mises 
yield condition yield condition 

0 0-1111 0-1067 

0-1 0-1086 0:104-2 

0-2 0-1011 0-0967 

0:3 0-0886 0-084-2 

0-4 0-0711 0:0667 

0-5 0-04.86 0-04-42 
0-6 0-0211 0-0167 
0-653 0-0044 0 

0-667 0 

0-697 7 0-0100 0-0135 
0-749 0-0276 0-0290 
0-814 0-0497 0-0482 
0-895 0-0783 0-0740 
0-962 -« 0-1028 0-0970 
1-00 0-1173 Oma 
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Résumé 


Dans ce mémoire nous étudions le probléme statique du choix des dimensions 
suivant: Comment doit varier l’épaisseur de point en point dans les plaques 
circulaires fléchies du type sandwich, pourque leur poids soit minimum, mais 

sans qu’il existe du fluage illimité? Comme hypothése de base, nous avons pris 
la condition de plasticité de von Miszs et la régle de fluage associée. Les plaques 
sont simplement posées ou encastrées le long de leur contour et portent une 
_charge répartie symétriquement par rapport a l’axe. Les résultats obtenus sont 
-comparés 4 ceux du probléme analogue, ou I’on utilise la condition de plasticite 
de Tresca. Nous avons présenté finalement quelques résultats numériques. 
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On Rank-Diminishing Operations and their Applications 
to the Solution of Linear Equations 


By EucEen EcrervAry, Budapest, Hungary’) 


Notations 


A=[a;;| | matrix composed of the elements a;; , 


Ax transpose of A, 

a column-vector, 

a* row-vector, 

<a,, ..., 4,> diagonal matrix, 
E unit matrix, 

0(A) rank of A, 

| A| determinant of A. 


Many special finite iterative methods for solving linear equations have been 
discussed by various writers?). The methods used apply for the most part to the 
special problems considered by each writer. 

In this paper we shall develop a general finite iterative method whiel 
consists in the repetition of a rank-diminishing operation and which will enable 
us to give a unified treatment for various iterative procedures occurring in the 
solution of linear equations. 

These iterative procedures can be classified in two categories: 

1, iteration of the rank-diminishing operation in order to obtain a dyadic 
decomposition of the matrix of a given system of equations, 

2. iteration of the rank-diminishing operation in order to reduce the whole 
vector space (the unit matrix) to the subspace (to the matrix of lower rank) 
which constitutes the solution of the given system. 

The first category includes 1. the triangular factorisation of GAuss-BANA- 
CHIEWIcz and its generalization, 2. the recursive method for the inversion of 
a matrix and its application to the solution of boundary value problems, 
3. the method of conjugate directions and the establishment of a general 
inversion formula including the inversion formula of STIEFEL-HESTENES. 

The procedures belonging to the second category may be interpreted as 
general projection methods of which the vector method of PURCELL-MOTZKIN 


1) Institut de Mathématique de l’Académie des Sciences de Hongrie. A posthumous paper, 
prepared for publication by P. Rézsa. 


*) See References [1]—[8] on page 386, 
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and the author’s method for solving a system of diophantine equations are 
particular cases. 


All these methods require the repetition of simple, elementary routines, 
hence they seem to be suitable for machine computation. 


1. Simple Rank-Diminishing Operations 


The well-known inequalities 
e(A) — 1 = 9A) — o(b c*) Se (A — bc*) < QA) + 0(b c*) = (A) + 1 


show that by the subtraction of a dyadic b c* from a given matrix A the 
rank of A increases or decreases just by one or it remains constant. It is easy 
to see that by this procedure the rank of A cannot decrease unless the left 
factor b is in the column space of A and the right factor c* is in the row space 
of A. But in the latter case b has the form A wu and c* has the form v* A, hence 
we have only to investigate the rank of A—A A uv* A, A being ascalar parameter. 
Let A = B C* be a basic factorisation (see e. g. [2]) of A. Then we have 


A—AAuv* A= B (E—AC* uv* B) C*; 0(E) = 0(A). 
Here B is invertable from the left, C* is invertable from the right and 
o (E—AC* uv* B) = 0(A) — 1. 
The value of the determinant 


|E—AC*uv* BJ) =1—Av*Au 
shows that 
o (E—AC* uv* B) = (A) — 1 


if and only if v* A u + 0 and A = (v* A u)-1. Hence we have the 


Lemma. By the subtraction of a dyadic b c* from a given matrix A, the rank 
of A decreases just by one if and only if b c* has the form 
Auv*taA (i) 
pA te 


where u and v* are arbitrary vectors subject only to the restriction v* Au + 0. 
Starting now from an arbitrary matrix A of rank 7, we can decompose it 
into a sum of r dyadics by the following iterative procedure 
A, Un 41 Yea A 
Ay Uys 


= (2) 


Ax = A, 


* 
taal 


where u,,, and v¥,, are arbitrary vectors subject only to the restrictions 
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vit, , Ay Up.1 + 0. Then the dyadic decomposition of A will be given by 


BLE (3) 


fete Ve u ve 
k+1 -k+1 
A= : ~ 


2. Multiple Rank-Diminishing Operations 


In this section we shall prove that the result A, of x successive simple 
rank-diminishing operations (1.2) can be obtained by the following multiple 
rank-diminishing operation 


A,=A—AU, (Vi AU,)1VEA, (1) 
where 
Ua 
U, = [tie tals VE | Ve ee 
vt 


Suppose that (1) is proved for x = 1, 2;..., &. Then we have by definition 


= = * ail cys: 
Ay 1 = Ay — Ay Up 41 Vii Ay Uys)? Vou Ay - 


Substituting A, from (1) this becomes 


A,.;=A-— AU, (Vi AU,) "VE A 


{A— AU, (VE AU) VE A} uy; v*,, {A — AU, (VE AU, VE A} 
Via {A AU, AVE A UVR TA) a 


‘Now we have to prove that the formula (1) for x =k +1 yields the same 
expression for A,,,. Writing U,,, and V#,, as partitioned matrices 


Les Ve AU, V¥ Au ‘ 

Nags = [U,, be 1 VE = | VE, AU : "see ss 
"hs Viv, AU, VE) A Uy 
the formula (1) for x = k + 1 takes the form ; 


V*¥ AU,, V* Au “17 V# A ; 
_A-[A te Amal]. é : aa y 4 (3) — 
v4.14 U,, DE ree 
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w 
“I 
Oo 


The inversion formula (2. 2) for partitioned matrices gives 


rea, VS Au. is a U,)-4: i 


Mi AU, ve AU, 0, 0 


—(V¥ AU) V* Au, . f 
et Akt. as [— vt, AU, (VE 4 U,)-4, 1] 


+ —_,—_,— Oe FS 
Mra A Uy — MF,, AU, (VE AU) VE AU, 


and substituting this in (3) we see that both expressions for A,,, coincide. 
Thus we have proved the validity of (1) by induction. 
If the rank of A is 7 then A, vanishes and we have the following represen- 
tation of A 
v*¥ A 
A=[Au,,...,Au,] [vt Au) ro (4) 
v* A 


3. The Triangular Factorisation of Gauss-Banachiewicz and its 
Generalization 


~It is known [1] that any matrix can be decomposed into the product of 
two triangular matrices provided that the upperleft principal minors are non- 
singular. In order to obtain the triangular factorisation of ann < m matrix of 
rank ry by means of rank-diminishing operations one can proceed as follows. 
Choose in the general formula 


mi 
0 
u,=e,=|1 (*; vt¥=ft=[0...010...0]); ¢=1,2,...,7 (1) 
0 ae Cy m 
. 0 | 


Then A e, is the first column, f¥ A is the first row of A and fk Ae, = ay + 0 
_ by assumption. Applying (1.2) with these vectors, obviously we get a matrix 


of rank ry — 1 


F Ae, ft A 
“ isiAs 141 


Sf Ae, : 


the first column and first row of which vanish. 
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The first non-vanishing column of A,-is A, @g, its first non-vanishing row 
is ff A, and the scalar f¥ Ay €: = (441 422 — 2 Ay1)/44, + O by assumption. 
Hence we can form the matrix of rank r — 2 
A A, &, fp A, 
die fae oe 


which has two 0-columns and two 0-rows. 
Continuing in this way we find A, = 0 and the result can be written as 


follows: fe A 
S# Ay 
A a [A e,, A, e,, ches Say A, e,] <Ji A e,; ge ae oN or e,> 
ie Ae 


The triangular form of the factors is implied by the 0-columns and 0-rows of 
the matrices A,. 

_ This particular method of factorisation fails if any of the left-upper principal 
minors of A is singular. A more general treatment [2] is possible in which the © 
unit vectors e;, f* are taken in any order. In this case the final factor matrices 
will, in general, not be triangular, but will be triangular matrices with their — 
rows and columns in some way interchanged. 


4. The Recursive Method for the Inversion of a Matrix and its 
Application to the Solution of Boundary Value Problems 


The rank-diminishing operation used in the preceding section furnishes a 
particularly simple result when applied to the inverse of a non-singular matrix. 
We write the given non-singular m-th order matrix in the partitioned form 


ne a, 
ae = 
a" a, n 
and suppose | A,_,| + 0. Then we shall prove that 


A-} A; e, en ae Ate 0 
tee igi a catan ae a (1) 


Indeed, the elements of the left upper block are proportional to the second 
order determinants formed with the elements of A, ', hence in virtue of a 


poeorem of JAcoBI they turn out to be equal to the corresponding elements of 
Avs 


. 
. 
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The relation (1), in its original form, can be used to compute the inverse - 
of a minor matrix, when the inverse of the original matrix is known. Hence it 
can be applied e.g. to the solution-of a difference equation with modified 
boundary conditions if the original boundary value problem is solved [8]. 

On the other hand, the relation (1) can be written in such a form, in which 


_ it is used in the escalator method for the inversion of a matrix. Indeed we have 


mre 1 
A A, 235 4, A, 274, oly aan SP 0 1 
cm a2 n A- 
4 1 ifs n a (a, mye KU Ans a,,) , 
a Qin 1 1 
hence 
1 2 
A> = Ann a, es ef 
ae rs (a, aye A, 1 a,,) 
it 
Similarly 
* -1__ n —1 7 n —) -1 
€, A, = [— a AS, ’ 1} (a, nm a A,-y a,,) : 
Finally 
* —1 aa n —1 -1 
en ZN eo (Ann way ASS, a,,) : 


Substituting these expressions in (1) we obtain 


_—A- : 
cE rll My Ay. 0 wast "a [= a" Are y5 tl 
La = ~ (2) 


a—arA'.a 
a a OFo-0 n-1 on 


nn 
n 


and this coincides with the well-known formula of the escalator method [3]. 
Thus it turns out that the escalator method furnishes the triangular factorisation 


~ of the inverse. 


5. The Method of Conjugate Directions. The Generalized Inverse. 
The Inversion Formula of Stiefel-Hestenes 


The dyadic decomposition (1.3) assumes a particularly simple form if one 
uses two sets of vectors w,,..., u, and v¥,..., vF which are conjugate with 


respect to A, 1.e. 


vt Au,=0 for. 1+7; vf Au, +0. (1) 
: In this case the matrix [v¥ A u,] in (2.4) will be diagonal and the formula 
~ reduces to anes: 
r,Au,vrA 
=) — Z 
A zen ve Au, - (2) 


” 


hoo 


g 


aA Yee ay 


This decomposition can be used for the computation of a generalized 
inverse of an arbitrary (rectangular) matrix A of rank 7. A generalized inverse 
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is defined as a common solution of the equations [5] 
AXA=A and XAX=X. (3) 


Now it is easy to see that in consequence of equations (1) the matrix of rank 7 


* 
Uy, Vz 


IE = ACV SD eA ee (4) 


1 


satisfies equations (3), thus it is a generalized inverse of A. 
If A is an n-th order, non-singular matrix, then multiplication of both sides 
of (2) by A-} from the left and from the right gives immediately 


* 
u, VE 


Wn 
At= ) —___— 
Dy vk Au, 
and this is the basic inversion formula of STIEFEL-HESTENEs [4]. 


6. Geometric Interpretation of the Rank-Diminishing Operation and 
Simultaneous Solution of two Systems of Equations 


The resulting matrix A — A u (v* A u)—! y* A of a rank-diminishing opera- 


tion can be factorised as follows 


Au v* 
ee ve 


an) A or A(E- era) 


v® Au 


Evidently E — (v* A u)-! A u v* is a projector, hence the factorisation on the 
left admits the following geometric interpretation : 


By the rank-diminishing operation (1.1) all column vectors of A are projected 
parallel to A u on the hyperplane v* « = 0. 


A similar argument shows that 


by the same rank-diminishing operation (1.1) all row vectors of A are projected 1 


parallel to v* A on the hyperplane x* u = 0. 


This observation gives us the tool for solving simultaneously two systems 


of linear homogeneous equations by a projection method (provided that some 
restrictive conditions are satisfied). 


Consider the following two homogeneous systems 


a* X=0 Xb, = 


a ae one 


oN a a hha 8 i 
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es) 


where a¥, ..., a," are linearly independent »-dimensional row-vectors, b,,*..,b, 
are linearly independent n-dimensional column vectors and X is the unknown 
n-th order matrix. If we succeed in finding a solution matrix of rank n —r then 
the columns of X furnish a complete solution of the system (la) and the rows 
of X supply a complete solution of the system (1b). 

We can find a matrix X, of rank n — 1 satisfying the first pair of the 
equations (1) by the following rank-diminishing operation 


Dg apne Bae 


at b, 


provided that af b, + 0. We have obviously a¥ X = 0; X b, = 0. 
Suppose now that a¥ X, b, + 0. Then the matrix 


X, b, af X, 
at X, by 


has the rank m — 1 and satisfies the first two pairs of the equations (1). 
Continuing in this way and supposing that the scalars a*,, X, 6,4, 
(k = 0,1,...,7— 1) are different from zero we arrive finally at a matrix 


X41 b, ay X41 


* 
ax X,_, 5, 


Xx, oa pe 


which has the rank m — ry and satisfies both systems (1). Thus its columns, 
or its rows, constitute a complete solution of the two systems. 

These results will be applied in the following sections to obtain the solution 
matrix of one of the given systems in a computationally convenient form. 


7. Solution of an Arbitrary Homogeneous System by a Projection 
Method. The Vector Method of Purcell-Motzkin 


Two subspaces of the n-dimensional vector space are given: the subspace YU 
spanned by the linearly independent column vectors of the matrix 


A= |a,, ...., G,] (1) 


and the subspace 8, defined as the intersection of the hyperplanes 
bt he 0, BPH 0, anny OF =P (2) 


(the b#’s are not necessarily linearly independent). 

It is necessary to find a basis for the intersection of W& and 8. 

Let b* be the first vector in the sequence (2) for which bX A + 0 and 
choose a non-vanishing element b* Ae, of b} A. Then by the rank-dimin- 
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Ae, by A 
My 
cr ae. 


My 


we obtain a matrix A, which has the rank 7 — 1 and satisfies the equations 


bF A, = 05) 9S 12 es 


ishing operation Are 


Furthermore it satisfies A,e, = 0, hence its “,-th column vanishes, con- 
sequently its remaining 7 — 1 columns are linearly independent. 

Let b* be the first vector in the sequence (2) for which b% A, + 0 and 
choose a non-vanishing element bx A, e,, of b* A,. The matrix 


ae an 


He Ve 


b* A, e 
Ve He 


A, = A, 


has the rank v — 2, it satisfies the equations 


bt A,=0, v=1,2,... 5%. 


Furthermore it satisfies the equations A, e,, = 0, A, e,,, = 0, hence its w,-th and 
{4.-th columns vanish, the remaining 7 — 2 columns being linearly independent. 

Continuing in this way, we arrive finally at a matrix of rank v — o (or to the 
0-matrix if the intersection of {{ and % is 0-dimensional) 


* 
A Pm: A A,_ 1 ing by A o-1 
ca o-l bx A e 
Ae 1 “Uo 


which satisfies all the equations 
yee: Wiad Pairs 


The w,-th, m-th, ... , 4,-th columns of A, vanish and the remaining 7 — o 
columns constitute a basis for the intersection of 2 and B. 
The method described here admits various applications. First, it stones an 


algorithm for ‘purifying’ a given set [a,,..., a,] of vectors from components — 


lying in the directions b¥, ... , b*. 
Furthermore if Q is the whole vector space, i.e. if A is a non-singular 


matrix (most conveniently the unit matrix) then this algorithm furnishes an 
independent set of m — o solution vectors for the general homogeneous system | 


(2). Moreover the algorithm eliminates automatically such equations (2) which 
are consequences of the preceding ones. 
In the special case when 


bt 


en 


pm Rt § as 


PR A a ao 8 me cee 


— 


- 


=A 
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is an mx m-+1 matrix belonging to a regular inhomogeneous system 
(written in homogeneous form) and we choose A to be the m + 1-th order unit 
matrix then our algorithm reduces to the ‘vector method’ of PURCELL and 
MOTZKIN [6]. 


8. The Solution of a Homogeneous System of Diophantine Equations [7] 


Let now 4, a, ... represent n-dimensional vectors the components of which 
are integers. We shall denote by A(a) the greatest common divisor of the 
components of a. 

We consider the system 

a*¥ X = 0 


(1) 
a*¥ X=0, 


where X denotes the unknown m-th order matrix the elements of which are 

integers, and suppose provisionally that a}, ... , a* are linearly independent. 
First we choose a vector b,, fixed once for all, such that a¥ b, = A(a#) and 

then compute the matrix of rank m — 1 

b, af 


* 
a} b, 


eee 


all the elements of which are integers. This matrix evidently satisfies the first 
equation a¥ X, = 0. 

a* being linearly independent of a, we have a¥ X, + 0 thus we can choose 
a fixed vector b,, such that a¥ X, b, = A(a¥ X,) and compute the matrix of 


rank n — 2 Wee iy ne Sas oe 
a 1 


af X;, 6, 


all the elements of which are also integers. X, satisfies af X, = 0 and a} X, = 0. 
Continuing in this way, we arrive finally at a matrix of rank n —7 


D. Gag sell Peat Chala Ge 
Rigi Mey ae x, i b, a a; Xia b, = A(aF X, -1) 
ap 


which satisfies all the equations (1). Any set of m —y columns, which are 
linearly independent, constitutes a complete set of solutions of the diophantine 
system (1). The general solution of (1) can be represented also in the form 


aa AL ey, 


where t denotes an n-dimensional vector whose components are arbitrary 


integers. 


ZAMP X1/25 
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If the vectors a*, ..., a* are not linearly independent then some of the 
vectors a¥ X,, k > 71 vanish. This means evidently that the corresponding 
equations a¥ X = 0 are consequences of the former ones and thus they must be 


neglected. 
REFERENCES 


[1] R. Zurmtur, Matrizen (Springer-Verlag, Berlin 1950), S. 248. 

[2] E. Ecervary, Uber die Faktorisation von Matrizen und thre Anwendung auf die 
Lésung von linearen Gleichungssystemen. Z. angew. Math. Mech. 35, 111 (1955). 

[3] A. Frazer, J. Duncan, and R. CoLiar, Elementary matrices (Cambridge Univ. 
Press 1938), p. 112. 

[4] M. R. Hestenes and E. StieFev, Method of conjugate gradients for solving linear 
systems, J. Res. Nat. Bur. Standards, 49, 409 (1952). 

[5] R. Penrose, A generalized inverse for matrices. Proc. Cambridge Phil. SOx Di, 
406 (1955). 


[6] E. Purcety, The vector method of solving simultaneous linear equations, J. Math. 


Phys. 32, 180 (1953). 
[7] E. Ecervdry, Auflésung eines homogenen linearen diophantischen Glerchungs- 
systems mit Hilfe von Projektormatrizen, Publ. math., Debrecen,4, 481 (1955-56). 
[8] E. EcervAry, Uber eine Methode zur numerischen Losung der Potssonschen 
Differenzengleichung fiir beliebige Gebiete, Acta Math. Acad. Sci. Hung. 77 (1960). 


Zusammenfassung 


Es sind in der Literatur (siehe [1] bis [8]) verschiedene iterative Methoden zur 


Lésung linearer Gleichungssysteme bekannt, die jedoch im allgemeinen nur bei den © 


speziellen Problemen angewandt wurden, welche die einzelnen Verfasser betrach- 
teten. 

In der vorliegenden Arbeit wird eine allgemeine finite iterative Methode ent- 
wickelt, die aus der wiederholten Durchfiihrung einer rangvermindernden Operation 
besteht. Diese Methode erméglicht die einheitliche Betrachtung verschiedener Ver- 
fahren, welche bei der Lésung linearer Gleichungssysteme auftreten. 

Diese iterativen Verfahren kénnen in zwei Kategorien eingereiht werden: 
1. Iteration der rangvermindernden Operation zur Gewinnung der dyadischen Zer- 


legung der Koeffizientenmatrix eines gegebenen Gleichungssystems; 2. Iteration — 


der rangvermindernden Operation zur Reduktion des ganzen Vektorenraumes (der 
Einheitsmatrix auf den Unterraum, die Matrix niedrigeren Ranges), welcher die 
Lésung des gegebenen Gleichungssystems bildet. 

Die erste Kategorie enthalt: 


1, die Zerlegung der Koeffizientenmatrix in das Produkt von zwei Dreieck- — 


matrizen (GAusSs-BANACHIEWICcz) und die Verallgemeinerung dieser Faktorisation; 


2, die rekursive Methode fiir das Invertieren einer Matrix und ihre Anwendung — 


auf die Lésung gewisser Randwertaufgaben; : 

3. die Methode der konjugierten Richtungen und die Aufstellung einer allgemei- 
nen Inversionsformel, welche die Formel von St1EFEL-HESTENES enthalt. 

Die Verfahren der zweiten Kategorie kénnen als allgemeine Projektionsmetho- 
den interpretiert werden, welche die vektorielle Methode von PuRCELL-MOTZKIN 


und die Methode des Verfassers fiir die Lésung diophantischer Gleichungssysteme 
als Sonderfalle enthalten. 


(Received: February 8, 1960.) 
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§ 1. L’évaluation par défaut de Payne-Weinberger 
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Nous étudions la premiére valeur propre A, d’une membrane vibrante 
‘couvrant un domaine G du plan et élastiquement liée le long de son contour I’. 


(La fréquence fondamentale est VA.) La fonction propre fondamentale u(x, y) 


satisfait A Aw + A, u=0 et wu > 0 dans G, du/on + R(s) u = 0 sur I’ (s est la 


longueur de l’arc, n la normale extérieure). - Aux points de J” oti A(s) = co 


I 


nous avons la condition « = 0 (contour fixé); la ot A(s) s’annule, le contour 


est dit libre. 


Une borne supérieure pour A, peut étre obtenue, a partir de toute fonction 
u(x, y) continue et suffisamment réguliére dans G, grace au principe de Rayleigh 


v) + f &(s) v? ds 


= } 
Pe dA 


4 = Min, 


1) Institut Battelle, Genéve, et Ecole Polytechnique Fédérale, Zurich. 
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ou D(v) désigne l'intégrale de DIRICHLET iG + v;) dA et dA = dx dy Vélé- 
G 


ment d’aire. Le minimum est atteint par la premiére fonction propre u(x, y). 

L. E. Payne et H. F. WEINBERGER [6]2) ont obtenu des évaluations par 
défaut pour 2, en comparant le probléme donné a celui de cordes vibrantes 
homogénes. Cette idée d’utiliser des problémes auxtliaires a une dimension est ala 
base du présent travail; comme on le verra au § 9, elle est intimement hée au 
principe de Thomson (problémes de la capacité et de la rigidité a la torsion, 
cf. [8]). 

Pour simplifier l’exposé, je résumerai la méthode de PAYNE-WEINBERGER 
dans le cas d’une membrane fixée: k = oo, c’est-a-dire wu = 0 sur J’. — Sur tout 
segment horizontal y,,: y = Yo, @ < x < 6 dans G, de longueur 6,,= 6 — a et 
a extrémités sur J’, ils considérent le probléme de la corde vibrante: 
t"(x) + v f(x) = 0, f(a) = f(b) = 0; », = 27/0), étant la premiére valeur propre 
de ce probléme auxiliaire, on a le principe de minimum 


Une inégalité analogue est valable pour les segments verticaux dans G, a 
extrémités sur J’. - En sommant toutes ces inégalités, on obtient 


Eile dA = Di) So? (rs + 7) [[weaa, 


ou L, est la borne supérieure des longueurs des segments horizontaux dans 
G, L, celle des segments verticaux; d’ou 


A, 2a? (Fe + sa (1) 


/ 

r ‘ ? : 
PAYNE et WEINBERGER traitent, d’une facon analogue, la membrane elas-— 
tiquement liée. | 


*) Les chiffres entre crochets renvoient a la Bibliographie, page 412 


Py Sen ee 
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§ 2. Méthode plus générale: comparaison a des cordes vibrantes 
inhomogénes 


Le raisonnement ci-dessus peut étre modifié en considérant, sur 24. Ut 
probléme auxiliaire 


P(x) + va(x, vo) (x)= 9, f(a) = f(b) = 0, 


a(x, y) étant une fonction choisie de x et y; le principe de minimum pour la 
premiére valeur propre 7,(Vo) est 


¥%4(¥0) = Minyy —-*—————. avec’ o(a) = v(8) = 0; (2) 


d’ou 
b b 


” ” 


| ile, Yo) de 2 ry(y0) | x(x, Yo) w(x, Yo) dx 


a 


a 


De méme, choisissons une fonction A(x, y) et considérons, sur un segment 


vertical y,: x= %, a< y <8, a extrémités sur I’, le probléme auxiliaire 
ge" (y) + o B(x, v) gly) = 9, g(a) = g(b) = 0; le principe de minimum pour 0,(%9) 
est 


foro 
Mi : 
0,(%9) = Min, Bigs oe 
Ja Xo, V 
u , ‘ 
dot 
b 
[W(x ») 4 1(% oo) Xo, V) U(X, y) ay « 
a 
En sommant toutes ces inégalités, nous obtenons 


Bate ls ee 


wey 


ly 


an 


FETT ROK 


MTs D(w) = [] rly) ale, ») + orto) Ble, v)] waa, 


G 


wis 
° 
fava 


[[ x — 00) a€e, 9) — onl) Ble, vl wt ad 2 0 
G 
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par conséquent, nous avons la borne inférieure 
dy = inf, peg [rrl) a(x, ») + o4() BOE, »)1- (3) 


Remarque. Si une horizontale y = Vy a, dans G, plusieurs segments distincts 
Yy.1 et Yy,2 le lecteur rectifiera de facon évidente la simplification excessive de 
notre notation »,(y): on aura une valeur de »(yo) sur y,,;, une autre sur y,,9- 

En spécialisant «(x, y) = B(x, y) = 1, on obtient une évaluation légérement 
meilleure que (1): 


1 1 1 1 ‘ 
A, 2m? inf, wee (se - a PE (ar ~ +z) : (1’) 


§ 3. Une expression du principe de maximum 


L’inégalité (3) est d’un usage malaisé: « et # étant choisis, on ne peut pas 
en général déterminer exactement les valeurs propres », et o, des problémes 
auxiliaires 4 une dimension. 


3.1. Retournons donc les problémes auxiliaires untdimensionnels! 


Choisissons, pour chaque y,, une fonction continue et deux fois dérivable 
f(x) = f(x, Yo) > 0; puis déterminons «(x, y 9) telle que f(x, vp) soit l’extrémale 
du probléme de variation (2): nous devons donc prendre 

fx(%, Vo) . 


V4 (Yo) &(%, Yo) = — “f(% Vo) 


il est bien entendu que la fonction auxiliaire f(x, y) n’est pas nécessairement 
continue en y. Choisissons de méme une fonction g(x, y) > 0 continue en y 
(mais non pas nécessaivement en x!) et deux fois partiellement dérivable par 
rapport a y; et prenons 0,(x) B(x, y) = — g,,/g; nous avons alors, au lieu de (3), — 
Vinégalité 


ji int (2 — fur). (4) 


3.2. Démonstration générale pour une membrane élastiquement liée 


Choisissons [3] dans G deux fonctions positives: 
f(x, ¥) continue en x et deux fois dérivable par rapport a x, g(x, y) continue en y 
et deux fois dérivable par rapport a y, satisfaisant sur I” a la condition 


fon gp: On = MS) (5), 


ou n est la normale extérieure, et le coefficient k(s) a été défini au § 1. ; 
i 
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On a alors la borne inférieure (4). 


Liberté essentielle (en vue des applications) : f peut étre discontinue en y, et 
g peut étre discontinue en x. 


Démonstration. Comme f + 0 dans G, ona 


et de méme 


g, og, 2 
uy + ae = i 1?) + (u, tes u) 
v 


d’ou, en ens iaae 


en tenant compte des égalités 


~~ - ) ) 
)+ pH) ) wds = Ay || u2dA , OY a Peo CLONAL easy 
nous obtenons 


i (4. + fee + Sur) wd = Diu) + f Rs) w® ds + I ee + $v) wdd 
E 4 é 


= a+ (fe 3 tae (le fe tan v0, 
G 


d’ou (4). 
Si l’on a choisi f(x, vy) = B(y) u(x, y) et g(x, vy) = P(x) u(x, y) (les fonctions 


Pet Y étant + 0 quelconques, continues ou non), le membre de droite dans 


(4) devient — Au/u = A,, on a donc l’égalité; nous avons donc le principe de 


maximum 


sia \ df al 


Ae ST PIG 


- Tove fu) 
A, = MaXnoix de fetg inf,,, y)EeG (- f g C 
satisfaisant (5) 


(7) 


3. 3. Remarques 


3.3.1. La borne inférieure fournie par (4) sera particuliérement bonne si l’on 
choisit f et g en sorte que (— f,,/f — &,,/g) soit constant (ou presque constant) 
_ dans G. 


392 JosepH HERSCH ZAMP 


3.3.2. Contrairement A ce qui se passe pour le principe de RAYLEIGH, la 
condition (5) est d’autant plus restrictive que le contour de la membrane est 
plus libre. Il devait bien en étre ainsi: le Maximum diminue. Si le contour est 
fixé (w = 0 sur J’), on a k = ov et la condition disparait ; on choisira alors avec 
avantage, sur I’, f = g = 0. 


3.3.3. Sil’on choisit f(x, y) = g(x, y) > 0 (donc fonction continue de x et de 
y), la condition (5) devient df/0n + k(s)f 20 et Vinégalité (4) donne 
A, = infg (— Af/f); c’est essentiellement l’une des bornes de BARTA-POLYA 
pour A,: 

Si f > 0 dans G et Of/On + R(s) f =O sur I, ona 


inf,(- a is An sup;(—“"). 


Démonstration directe: 


ap 


[fu (at+apaa = |] wAt— fu) dA = o (u So — 7 SE) as 
= f Hs) (—uf+fu)ds=0, | 


d’ot les deux inégalités ci-dessus. Ainsi, toute fonction réguliére /, positive — 
dans G et s’annulant sur J’, fournit des bornes inférieure et supérieure pour la 


premiére valeur propre de la membrane a contour fixé (w= 0 sur I’).—Ce — 
dernier théoréme est cependant d’une application pratique malaisée. 


3.3.4. Les fonctions », «(x, y) = — j,,/f et o, B(x, vy) = — g,,/g peuvent par- 
faitement étre par endroits négatives. Les cordes vibrantes auxiliaires que nous 
considérons (cf. §2) peuvent donc avoir par endroits une masse spécifique 
négative! 


§ 4. L’expression vectorielle du principe. - 
Coordonnées cartésiennes et coordonnées polaires 
4.1. Passage de l’expression précédente a la forme vectorielle 
A partir des fonctions f(x, y) et g(x, y) choisies au § 3.2, définissons le champ 


vectoriel 
Ad he tiie by). : 


il satisfait a la condition sur le contour 
pon < Rs), (5! 


ou de nouveau nm est le vecteur unitaire indiquant la normale extérieure. — 


| 
: 
| 
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On calcule sans peine div p — p? = — fad Bile done 
A, = inf, (div p — p?) . (4’) 


Liberté importante pour les applications: de champ p n'est pas nécessairement 
continu, sa premiére composante f, doit étre continue en x et dérivable par 
rapport a x; le seconde /, doit étre continue en y et dérivable par rapport a y. 

La condition / > 0, g > 0 dans G se traduit ici par | p | < oo. 

Ce passage a la forme vectorielle correspond au passage classique d’une 
équation de StuRM-LIOUVILLE a une équation de RiccaTI. 


4.2. Démonstration directe 
0 S (grad w+ u p)? = grad? u + p: grad (uw?) + 2 p? 
= grad? u + div (vu? p) + wv (— divp+ p’), 
avec égalité a gauche si et seulement si p = — grad u/u; d’ou par intégration 


0<D(u)+ f uw p-nds+ |] (—divp+ p) wad 


< D(u) + f k(s) 2 ds + |] (— div p + p) wd = ff ts — divp +p) wd; 


d’ou (4’). 
Si lon a choisi p = — grad w/w, la condition (5’) est satisfaite et le membre 
de droite dans (4’) devient —Au/u = A, , on a donc l’égalité; d’ot le principe 


_ de maximum 


: 3 2 
| A, = MaXgpoix de p inf, y)EG (div p 7 Pp ) : | (7’) 
satisfaisant (5’) 


4.3. Remarques 


4.3.1. Pour obtenir une bonne borne inférieure (4’), on construira de préfé- 


rence un champ p réalisant div p — p? = const dans G, et p- n = k(s) sur I’ 


a 


er ee N 


(il en est ainsi pour le champ extrémal — grad w/w) ; si k = oo (membrane a con- 
tour fixé), la condition (5’) tombe; mais on fera un usage judicieux de sa 
liberté en choisissant p singulier sur 1”. 


4.3.2. Les discontinuités admises pour le champ vectoriel p sont les mémes 
_ qu’admet le principe de Thomson (électrostatique, rigidité a la torsion, etc.). 
_ Pour les applications, la régle suivante suffira: lors du passage a travers une 
~ courbe analytique (dans G), la composante normale de p doit étre continue, la 
_ composante tangentielle peut étre discontinue. 
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4.3.3. L’expression vectorielle du principe n’est pas restreinte aux coordon- 
nées cartésiennes. En outre, elle permet d’abréger légérement les calculs. En 
revanche, l’énoncé donné au § 3 est plus intuitif et proche d’une interprétation 
physique (cf. § 10); il laisse aussi mieux apparaitre les relations qui existent 
avec la méthode de PAYNE-WEINBERGER [6] et les inégalités de BARTA-POLYA 
($:3.3:3): 

4.4. Coordonnées cartésiennes 


p={br bo}, Prilx, y) continue en x et dérivable en x, p2(x, y) continue en y et 
dérivable en y; satisfaisant 4 la condition p, 0x/0n + Pz dy/On < R(s). — Posons 


n) ¥ 
tales) = Fe — Bi 

op = Valores Aq <o inf, sec [xa(%; y) + 49(*X, yy). (8) 
Yolx, y) == aR LE 5 


Si l'on choisit 7,(x, Vo) et Y2(%o,y), ce sont 1a des équations du type de 
RiccaTl pour p,(x, Yo) et £2(%, y). Nous avons ici rejoint Vénoncé de M. H. 
PROTTER [9], mais avec des conditions de continuité plus générales (cf. § 5 cl- 


dessous). La ot 0u/x = 0, on choisira avec avantage p, = 0; 1a ot Ou/dy = 0, 


on choisira p, = 0. 
Champs particuliers (en vue des applications) : 


Si pi(%, Yo) = ¢, alors 471(%, Vo) =—.¢*; 
Si P1(%, Yo) = & tg(h x), alors 4%, Yo) = RP; 
si P(x, Yo) = — A ctg(k x), alors 44(%, Yo) = k2; 
Si £,(%, Yo) = — 2 tgh(x x), alors ,(%, Yo) = — 235 
si p1(*, Yop) = — xctgh(x% x), alors 4,(%, Yo) = —-x?; 
cas intermédiaire k = x = 0: 
, 1 
si p,(%, Yo) = — pit alors: 71(%, V5) =0. 


4.5. Coordonnées polaires 


P = {hi bo}, ict py = composante radiale, p, = composante azimutale. 
p,(7, 9) continue en 7 et dérivable par rapport a 7; 
(7, 0) continue en 6 et dérivable par rapport a 0; 

satisfaisant a la condition p, Or/On + 1 p, 00/0n < k(s); définissons 


walt, 6) = + 5 (bs) — 83 


aa : ; alors A, 2 inf, » <¢ [m(7, 8) + xol(7, 9]; (8 
x2(7, 0) = Fie RHE 
on a l’égalité si et seulement si p,=—u,/u et py = — 1/r w/w. 
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Champs particuliers 


Si, Z,(r) étant une fonction cylindrique d’ordre q; 


1 A 2 
Py(r, Oo) = — Fike : Z(k 7), alors’ y,(r, 0.) = h2 — 
Q 


dy v2 , 
si p,(%, 6) = m tg (q 0), alors: .¥4(%o, @) = cake 
0 
si Pal’ 8) =— 4 ctg(g 6), alors alr, 6) = 
0 


§ 5. Un travail de M. H. Protter 


Dans un «Technical Report» [9], non encore publié et dont j’ai eu connais- 
sance récemment, M. H. ProtTrerR a énoncé et démontré en 1958 le théoréme 
suivant pour la membrane fixée: 

«Sotent P(x, y) et Q(x, y) deux fonctions continiiment différentiables dans le 
domaine G. Alors, tout nombre A tel que P.+ Q, — A> PP + + Q? dans tout G, est 
une borne inférieure pour A,: 4, = A.» M. H. PRotter indique aussi la générali- 
sation a une membrane élastiquement liée. 

En d’autres termes, A, = inf, (P, + Q, — P? — Q?), ce qui n’est autre que 
la formulation de notre principe de maximum donnée au § 4.4 (poser #, = P, 


Pp. = Q); a cela prés, cependant: comme nous l’avons vu, les conditions de 


continuité (et de dérivabilité) imposées par PROTTER a P(x, y) et Q(x, y) sont 
plus strictes qu’il n’est nécessaire: cela restreint les applications du principe, et 
c’est pourquoi les exemples qu’il indique (cf. aussi [5]) ne sont pas trés concluants 
(la meilleure méthode connue a comparer, pour le cas du losange, a celle qui 
nous intéresse, est certainement la convexité de P6LYA-SCHIFFER [7]). I] n’en 
reste pas moins qu’a M. H. Protrer revient le mérite d’avowr découvert, le premier, 
une inégalité du type que nous étudions ict; il y est parvenu bien que, de toute 
évidence, il n’ait pas eu connaissance du travail, légérement antérieur, de 
L. E. Payne et H. F. WEINBERGER [6]. Celui-ci ouvrait la voie et c’est vrai- 
ment, comme le montrent les §§ 1, 2, 3 et 4, le point de départ du présent 


- mémoire. 


§ 6. Méthodes analogues dans le plan 


Au lieu de choisir deux directions privilégiées (orthogonales) « et yy, comme 


~ nous l’avons fait aux §§ 2 et 3, on peut également en choisir davantage. 


A 


Aa ie 


ee 


6.1. Par exemple, considérons ¢rois directions privilégiées §, 7, ¢ formant 
_ deux a deux des angles de 120°. Je désigne de nouveau par u(x, y) la fonction 
_ propre fondamentale de la membrane G élastiquement liée le long de son 
Be sctour I’ Au +4, u=Oet u > 0 dans G, dujdn + R(s) w= 0 sur I. 
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Notation: Pour toute fonction f dans le plan, le symbole /, désignera la déri- 
vée partielle de f dans la direction §. On vérifie facilement que uz+ u, + uE= 
3/2 grad? w. 

Considérons dans G trois fonctions Positives: 

{ = fonction continue en & et deux fois dérivable par rapport a & le long de 
toute paralléle 4 &; g = fonction continue en y, etc.; h = fonction continue en 
f, etc.; satisfaisant en outre a la condition 

he (OF bn On h- of a: 


R(s) « (5") 


3 
“f = on ee on i. on 2 


En chaque point de G, nous avons trois relations du type (6), cf. § 3.2; en 
les sommant et intégrant, on obtient 


aes Ie San Ae 2 = (0: 
a 
donc, sous la condition (5"), 
Tiaras fee g, is 
pes ss mm c " 
A, 2 = inf, ( ; - ale (4) 


On a l’égalité en particulier si f= g=h= u. 


6.2. Cas limite: toutes les directions sont «privilégiées», a poids égaux. Pour 
sumplifier, restreignons-nous aux membranes liées: k = co, u=0 sur I. 

Spécifions une direction § par l’angle 6 (0 < 6 < a); pour chaque valeur — 
de 6, nous choisissons une fonction /®, qui doit étre, sur chaque segment de 
direction @, continue et deux fois dérivable. Alors 


2 ee Tt oe 
Ay = ip Bc, (<a 5 am) ag. (4"") 


La démonstration est analogue a la précédente et tient compte de l’égalité 


‘ 
MI 

7 
; 
, 
‘ 


Tu 


"{ OW V2 
i: Car do = 5 grad? uw . 
0 2) 

\ 


En particulier: Appelons L(P, 6) la longueur (finie ou infinie) du segment 
rectiligne dans G, passant par un point P&G et de direction 6. Si ce segment | 


; 


pals oaks ee WA A+ = = 


= 


vr 
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est fini, on peut y choisir / = arc de sinusoide, en sorte que 
(6) 
fx) (0) or 
fo _ LE’ 
d’ot 
Fee (* ao 
4 = 2 i pee | : 
220 trea | LP, 0) ° (9) 


0 


De facon encore plus grossiére: soit L(6) la borne supérieure des longueurs 
des segments dans G et de direction 0; alors 


ra @ cag 

Gy Pe ' 

Pelee co) 
0 

Si G est une bande infinie de largeur a, L(0) = a/sin 0, d’ot 


; 2 
Ay = 2a as 


ce qui est la valeur exacte. 


§ 7. Expressions du principe pour l’espace a trois dimensions 


Nous considérons le probléme classique caractérisé par les équations 
Au + A, u=0 et u > 0 dans G, du/dn + k(Q) w= 0 sur I’, ot 


gh P haart aS 
= oxt Taye T at 


et k(Q) est une fonction donnée sur la frontiére I’ de G. n est de nouveau la 
normale extérieure. 


7.1. Transposons notre principe de maximum tel qu’il est énoncé au § 3.2. 
Choisissons dans G trois fonctions positives 


f(x, y, 2) continue en x et deux fois dérivable par rapport a x , 
g(x, y, 2) continue en y et deux fois dérivable par rapport a y , 
h(x, y, 2) continue en z et deux fois dérivable par rapport a 2, 


satisfaisant en outre, sur J’, a la condition 


On a alors la borne inférieure 


h 
Ay = inf (— ea (4%) 
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La démonstration utilise l’égalité D(u) + fp R(Q) w dS = A, I wdaV, ou 
r G 


dS désigne l’élément de surface, dV l’élément de volume et D(u) Vintégrale de 


DIRICHLET Hh grad? uw dV. Comme au § 3.2, on obtient en intégrant la relation 
iG 


(6) et en lui ajoutant les deux relations analogues en g et y, fi .ev ze 


ii] (4+ ies ans ee a ss) ydV 20, 
vs 


d’ot (41). 
On a l’égalité si f(x, y, z) = Bly, z) u(x, y, 2), g(%, ¥, 2) = Wi (2, x) a ey 2) et 
h(x, 9, 2) = 2(«, y) W(x, 9, 2). 


x 


7.2. Au lieu de comparer notre probléme a trois dimensions a celui de 
cordes vibrantes (inhomogénes), orientées dans 3 directions, nous pouvons 
également le comparer a celui de cordes vibrantes (inhomogénes) ortentées dans 
une direction et de membranes vibrantes (inhomogeénes aussi) situées dans les plans 
perpendiculaires a cette direction. On rejoint alors la spécialisation suivante du 
Shes 

Prenons z comme direction privilégiée; posons alors, dans le § 7.1 ci-dessus, ~ 
f(x, Vv, 2) = g(x, y, 2); nous avons a présent deux fonctions positives: 

f(x, ¥, 2) continue en x et y et deux fois dérivable par rapport a x et y, — 
h(x, y, 2) continue en z et deux fois dérivable par rapport a z, satisfaisant, Thy 7 
étant, dans un plan horizontal, la normale extérieure a la section plane de I’, 
la condition 


7 Og, ara eaaietaea BY) 


sous cette condition, 


Ay Bink, (— eet fw _ Ae), 4%) 


On a l’égalité si f(x, y, z) = O(z) u(x, y, 2) et h(x, y, 2) = Q(x, y) u(x, y, 2). 


7.3. Il y a évidemment encore d’autres possibilités pour construire des 
problémes auxiliaires 4 une et A deux dimensions (par exemple de facgon 
analogue au § 6). Nous allons encore considérer la comparaison du probléme — 
donné a celui de membranes vibrantes inhomogénes et respectivement paralléles aux 
trois plans (x y), (y 2), (zx). 

Considérons trois fonctions positives dans G: 


TX, O82) continue en x et y et 2 fois dérivable par rapport a x et y, 
g(x, ¥, 2) continue enyetzet 2 fois dérivable par rapport ay et z, 
h(x, y, 2) continue en z et x et 2 fois dérivable par rapport a z et x, 
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satisfaisant sur J’ a la condition 


ek) eee te OK on 
iv rs al 9 a ml: — oe VI 
fo OMyy On g on,, —On h Om, on B40) (5) 
alors 
A = te inf, (— fat toy = Suy + S22 re nest Peo) (4¥4) 
4 si g h : 


Démonstration. En vertu de (6), 


fog ; hy u 


2 


us + ur + Petit. {22 ut ; ty tl = div, mae grad,, /) 
et deux expressions analogues pour g et /; d’ou en intégrant 


a 


[ff (22+ fot te 4 Bvt berg Bethe) yay 20, 
le 5 


dot (4%), 
On a Végalité si 
f(x, ¥, 2) = Ps(z) u(x, y, 2), g(%, v, 2) = V(x) u(x, y, 2) 
et 
h(x, y, 2) = Poly) u(x, y, 2) « 


7.4. Expression vectorielle. Les formules (5’) et (4’) du § 4 restent valables 
pour trois dimensions; ici, p = {p,, p>, ps}, py est continue en x et dérivable par 
rapport a x, p, est continue en y, etc. Mutatis mutandis, la démonstration du 
§ 4.2 reste valable, ainsi que les remarques du § 4.3. 


§ 8. Quelques exemples d’applications 


8.1. Théoréme. Soient G un domaine convexe du plan et o le rayon du plus 
grand cercle inscrit dans G. Alors 


4 S@A(G) S75, 


ou j) = 2,4048 est le premier zéro de la fonction de BEssEL d’ordre zéro, et 


x 


A, se rapporte a la membrane a contour fixé. Ces bornes ne peuvent pas étre 


_ améliorées: On a l’égalité, 4 gauche, si G est une bande infinie; a droite, si G 
e est un cercle. 


c 


Démonstration. L’inégalité de droite est une conséquence immédiate du 


- fait que J, est une fonctionnelle de domaine monotone. Pour démontrer l’iné- 
_ galité de gauche, cherchons ov sont les points de contact de la circonférence C, 
_ du plus grand cercle inscrit avec le contour I" de G. 


ee 
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Tout arc connexe « d’angle > a sur la circonférence C a contient au moins 
un point de J’. En effet, si tel n’était pas le cas, on pourrait translater le cercle 
d’une longueur finie 4 l’intérieur de G, puis on pourrait augmenter le rayon @ 
d’une longueur finie également, sans que le cercle quitte G; cela contredit la 
définition de @. 

Soit 6 l’angle du plus grand arc de C, sans point commun avec I’; nous 
venons de voir que f <a. Si B =a, G (étant convexe) est inclus dans une 
bande infinie de largeur 2 9; donc (monotonie!) A,(G) = (2/2 @)?. Si au contraire 
B <2a,G est inclus dans un triangle T circonscrit au cercle C,; A,(G) 2 A,(T) 
(monotonie) ; il reste seulement 4 montrer que A,(7) = (2/2 @)?. 

Pour cela, tracons les trois segments rectilignes (bissectrices) joignant le 
centre de C, aux trois sommets de T; on a découpé ainsi T = T, — T, — T;; 
désignons par ¢,, ty, t; les cOtés correspondants de T. Appelons &, la premiére 
valeur propre d’une membrane sur 7,, fixée le long de #; et libre le long des 
deux «coupures» qui bordent T,; A,(T) = min (&,, &, €,): car le plus petit &; 
n’est autre que la premiére valeur propre de la membrane T décowpée le long 
des trois segments intérieurs. (Cf. la «méthode de découpage» [2], [4]). Nous 
allons évaluer §; a l’aide du principe du § 3: prenons-y ¢; comme axe des %; 
choisissons dans T7;: f(x, y) =1; g(x, y) = sin (4/2 @); la condition (5) est 


satisfaite (sur ¢;, on a k = oo; sur les deux coupures, k= 0 et g > 0, g, > OF 


Oy/On > 0); d’ou en vertu de (4), &; = (7/2 @)?; donc 


. rs 2 
A,(G) 24,(T) 2 min (&,, §5, €3) 2 ace 
ce qui achéve la démonstration. 
8.2. Fréquence fondamentale d'une membrane en L (fig. 1), fixée le long de son 


contour. 
By ee 1 iat 


8.2.1. Pour calculer des bornes inférieures pour A,, nous appliquerons le 
principe de maximum sous la forme du § 4.4. Je partage le domaine G en trois 


régions I, II, III comme l’indique la figure, j’y définis Px(%, y) et Po(x, y) ainsi: 
Dans I: £,= —A, ctg( 

Pe = —h, ctg ( 
Dans I: p, = —he ctg(ky%), dod y, = ke; 


wal dot i= ry, 


R 
; Mt Xe = 2k. 
ky 4), dot ¥, = h?; : : 


I 


( 
Dans III: 4, -= cte(= m4), dott y, = (=); m\2 
P2 = — he ctg (Rp yo), Wot x2 = 3; ee () = 
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Pour que le champ p ainsi défini soit admissible, il faut que Pp, soit continue 
en x et , continue en y, d’ot la conditiun 


—h, ctg (Rk, a) = +k, ctg (Ry d). (10) 
20 7 - 
an? - 
18 = 
z f Ie 
16 " } 
| } 5 fas 
4% }—-§ 4 
a\ \ Neen 
=e } GZ + 
Py Ze 
12 20 
is 
| Lally se Ne ot 
<= 10 — ee ——— 
% st? 
85835 — = a me ---+-4-- 
0 02 fh 5 Op wen ae 10 
Figure 1 Figure 2 
Membrane en L, a contour fixé. La borne inférieure est fournie par la méthode proposée. 


x 


Conformément 4 une remarque faite au § 4.3.1, nous obtiendrons une 
meilleure borne inférieure pour A, en choisissant k, et k, en sorte que’ 
4x + %2 = const dans G, c’est-a-dire 

15 


2h = (=) 4 Re. . (11) 


a 
Les deux relations (10) et (11) déterminent k, et k, en fonction de a et 0. 
Lorsqu’elles sont satisfaites, on a la borne inférieure A,(G) = 2 kj. Posons 
encore 
: PERG, C= hy 


7 £ 2 pes ; tg &\2 
a EF» EE et ie) =2ee arr do tert Oe att) (#2). 
3 
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Je prends maintenant ¢ comme parametre; j’en tirerai d’une part S puis 
b/a; d’autre part, la borne a? Ap eee 


1 Cas a = ie 
ee 
igé.. i ee 
tg¢= -/2f—2 § * yen yore re 
BONS IS ES 
ie 
ee tee b n 
C=17, nréel; tg =ttghy; tghy= Vn? — 2 2& ree = 


8.2.2. Il est facile (méthode de Pérya) d’obtenir une borne supérieure pour | 
4,(G), donnée par la premiére valeur propre du rectangle de cétés a et a + 2 b; 
en effet, sa fonction propre fondamentale fournit, par transplantation, une 
fonction concurrente pour le principe de RAYLEIGH de G. Donec 

1 1 
eee | 
A(G) Sax (=. ae on 25 iE 

8.2.3. Cas limite bla > 0. Pour b=0, on sait que a? A, = 27; prenons 

&=z — 6; alors 


&, 
> 

VV 
bo 
a 

| 
2) 
2 


2 4 a eS 2nt— 42 af 


D’autre part, nous avons 


re 


aa, Sa [1 + (1+22)"] 208 — 4 a? 


Par conséquent: 


a h, Ad 2nt— 42 


La courbe donnant a? J, en fonction de b/a est inconnue; mais on connait sa 
tangente au point b/a = 0. 


8.2.4. Le rapport b/a varie de 0 a oo lorsque le paramétre € diminue de % 
a environ 2,07165 (valeur pour laquelle tgh "= 1 et n = oo). 
Pour € = z, nous obtenons bla = 0 et. a2 A, = 222: 
pour § = n/V2 = 2,22144: b/a = 0,5914 et x? < a? A, < 1,210 2°; 
pour € ~ 2,07165: bla = 00 et 8,5835 < aA, < n?. 
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La figure 2 donne les évaluations par défaut et par excés pour a? 4, en 
fonction de b/a. 


8.3. Membrane circulaire de rayon 1, fixée le long de sa circonférence et sur le 
segment réel -1 < x S 0 (fig. 3). 

Nous calculerons des bornes inférieures pour A, A l’aide du § 4.5. Partageons 
G en I et II par le cercle concentrique de rayon | @ |. Nous devons considérer 
séparément les cas 9 SO et 0 => 0. 

8.3.1. Bornes inférieures pour le caso < 0. 


Définissons les composantes radiale #, et azimutale p, du champ auxiliaire 
Pp ainsi: 


Dans I: ,(r, 6) = y+ Gp Jon) = +h ee dou y,(7, 6) = 
blr, 0) =0, ~d’ot x,(r, 0) = 0; -4 + y= Re. 
Dans II: £,(r, 6) = x - Z Zyo(Ro 7) = he {+ ; Es are |, 
dob m(r, 0) = 5 Zip =a Jy t+ B Nays: 


Pare F. 2 
Srja(*) = Vaz sinx; Ny (x) = V4 cos x; 
2 ; ie ae 
J-\2(*) Sed | Bees cosx; N_4(4) = =e sinx . 
é 


1 Ras 1 5 2 
Plt, 0) = — te dot XAT, 0) es Xi + Xe = Ry - 


Pour que ce champ p soit admissible, il faut que #,(7, 0) soit continue en 


E y = — 0, d’oti la condition: 
k I Z_112 (— Re a = i= hy Q) ; 
2| 2he@ Z1;2 (— he @) * Jo (— *; Q) 


ib ey eu. 


TOOT aoe 


a 
- 
g 


a 
a 


in oy oe cia 


Nous obtiendrons une meilleure borne (cf. § 4.3.1) en choisissant y; + 7% = 
const dans tout G, c’est-a-dire en prenant 


ky = hg; 
et en choisissant /, singulier en 7 = 1, c’est-a-dire 
Zi)2 (ky . 1) — 0) . 


Nous avons donc trois relations qui déterminent k,, k, et «: B en fonction de 


o. Posons 
—ko=v0, 
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nous obtenons 


ctg(k — v) = AG Bai 
Prenons maintenant v comme paramétre: a l’aide de (12), nous déterminons 


d’abord la borne k < VA, puis @ = — v/k. 
Lorsque v décroit de 7) = 2,4048 a zéro, 9 croit de — 1 a zéro et & croit de 


7) 4 az. Nous constatons que VA, = k pour 9 = — let pour 9 = 0. 


24 
=]0 -08 -06 -04 -02 0 02 0h 06 (08 S70 
Ci eae ae 
Figure 3 Figure 4 
(Cas @ > 0). La borne inférieure est fournie par la méthode proposée. © 


8.3.2. Bornes supérieures pour le cas 0 < 0. ; 


On sait [8] que VA, Sj,/7, ot” est le rayon conforme maximum. Ici, 9 < 0, 
le maximum ¢ est atteint en un point réel x (0 < « < a\2= hava représentation 
conforme fournit les relations entre o,xetr: 


GO ae ae et ee ee 
a? ~ Gay eee © "=" gag 


Prenons donc « comme paramétre, nous calculerons les valeurs correspon= 
dantes de 9, 7 et j,/7 > Vay 


D’autre part, il est aie (monotonie) que Aviecoy aie 0 = 7. 


a 
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8.3.3. Bornes inférieures pour le cas @ = 0. 


Définissons un champ auxiliaire p par ses composantes radiale p, et azimu- 
tale p, (cf. § 4.5). 


, 1 d 
pee ae) may ae Bes alt a te) 
dot yx, (7, = w= 


1 Ppt. if 5 
Pol’, By see he dot xalr, 8) =; i+ x%2.= ky. 


: —_ i ad 1 Z_x10(h 7) 
pe tir Ps. 8) Z112(Re 1) a Zyalke dee (+ 2a eA | 
, 1 
d’ot. y,(7, 0) = kp — ae 


Pour que ce champ p soit admissible, p,(7, 6) doit étre continue en 7 = @, 
d’ou la condition: 


hy | 1 Jo(hr a > f 1 Zapalheey \. 


hgo — Shl¥i 0) Pye  “Lilglta 0) 


‘Pour obtenir une meilleure borne, nous choisissons 7, + 72 = const dans 
tout le domaine, donc k, = hk, et f, singulier en 7 = 1, donc Z,;.(R2) = 0. En 
posant k @ = w, nous obtenons 


an 1 Jo(w) , 
ctg (k — w) ae ale (12’) 
Prenons w comme paramétre, nous pourrons calculer les valeurs correspon- 


dantes de la borne k = Vay et de 9 = v/k. 
Lorsque w croit de zéro a 7, = 3,8317, @ croit de zéro a 1 et k croit dew a 4}. 


On voit que VA, = k pour 9 = 0 et pour 9 = 1. 


8.3.4. Bornes supérieures pour le cas 9 = 0. 
A cause de la monotonie de A,, une borne supérieure est fournie par la 


; premiére valeur propre d’une membrane circulaire de rayon 1 — @, fixée le 
_ long de sa circonférence et sur un rayon: Vay <z/(1— 4). D’autre part, égale- 
- ment A cause de la monotonie, VA,(0) < VA,(1) = 7, = 3,8317. 


a 


“I 
c 


LOA oe led has ae 


8.3.5. La figure 4 indique les bornes inférieures et supérieures obtenues pour 


Vay en fonction de g. On remarquera que la borne inférieure (fournie par le 


, 


principe de maximum) se confond avec la valeur exacte pour 9g = — 1,9 =0 
et o= +1. 
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8.4. Comparaison de la fréquence fondamentale d’une membrane avec celles que 
l'on obtient en déformant linéavrement le domaine dans trois directions 


8.4.1. L’inégalité suivante est une conséquente immediate du § pe 

Soit G un domaine de l’espace a trois dimensions, de surface I’; appelons 
Jz la premiére valeur propre du probléme Au + 1 u = 0 dans G, w= Osun 
appelons /,, la premiére valeur propre de la membrane (fixée) couvrant l’inter- 
section de G par le plan z = 2, et A,, = inf,, A,,; définissons de fagon analogue 


ee et 1 pe alors ie ae dee OT As : (13) 


Cette inégalité est, trés simplement, l’une des transpositions possibles, pour 
trois dimensions, de ]’inégalité initiale de PAYNE-WEINBERGER [6]. On a l’éga- 
lité si G est un parallélépipéde rectangle, d’arétes paralléles aux axes Ox, Oy, O2. 


8.4.2. Soit G un domaine du plan, 4,(G) la premiére valeur propre de la 


membrane fixée qui le recouvre. Construisons, dans l’espace, le prisme infini G 
de section droite G. On sait que Az = 4,(G). Tout plan o non paralléle aux 


génératrices du prisme, coupe G suivant un domaine A, quin’est autre qu’une © 


déformation affine de G, orthogonale a l’intersection des deux plans. 
I suffit de choisir trois plans 0,, 62, 03, orthogonaux deux a deux, pour que 
Vinégalité (13) donne 


Ay(Ag,) + Ax(Ag,) + An(Ao,) S 2 An(@) « (13") 


2 


Cette inégalité est 4 rapprocher de la propriété de convexité de A, (cf. POLYA- 
SCHIFFER [7]). Si l’on choisit og paralléle aux génératrices du prisme, /,(A,.) est 
a remplacer par z?/L?, ot. L est la borne supérieure des longueurs des segments 
dans G et paralléles a la trace de o3. 


§. 9. Relations avec le principe de Thomson 


Pour montrer que le principe de maximum qui fait l’objet de ce travail est 
trés proche du principe de THomson, je vais rejoindre celui-ci, pour les pro- 
bléemes de la capacité électrostatique et de la rigidité a la torsion, en raisonnant 


5 


comme aux §§ 2-4: en comparant le probléme donné a des problémes auxiliaires 


a une seule dimension. 
9.1. Capacité électrostatique 


9.1.1. G est un condensateur dans l’espace, J’, et I’; ses électrodes (surfaces 


' 


supposées suffisamment réguliéres). Le probleme de DrricHLet: Au = 0 dans Gq 


u = Osur I%, w= 1 sur I, permet de définir la capacité C par 


40C = Diu) = |] grad* wav = eee 
G ES 


' 
4 
, 


-- 


- 


ol at ain a oat: se ae 
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9.1.2. Construisons des problémes auxiliaires 4 une dimension: choisissons 
dans G trois fonctions «(x, y, z), B(x, y, 2), y(x, ¥, 2). Sur chaque segment f || O x 
dans G, a extrémités sur I, stpposons résolu le probléme aux limites: 
I(%, Yo, Zo) = 0 aux extrémités, f,,(%, Vo, 2) = a(X, Vo, Zo); la solution f(x, yo, 2) 
de ce probléme est aussi l’extrémale du probléme de variation: 0(%, Vo, 29) = 0 


. 


aux extrémités, | (vt + 2av) dx = min. La solution uw du probléme initial 


x 4 
dans l’espace est concurrente, sur t, pour ce probléme auxiliaire ; donc 


[ (uh + 2am) dx =/(f 4+ 2af) dx. (14) 


ny 


Tv 


Si t était a extrémités sur J\, on aurait imposé f= 1 aux extrémités; si + 
reliait J) 4 J, on aurait imposé f = 0 a l’extrémité sur I, f = 1 a l’extrémité 
sur I}. 

Appelons g(x, V, 29) les solutions des problémes analogues, a une dimension, 
sur les segments || Oy, avec g,, = B(%o, V, 2); et (Xo, Yo, 2) les solutions des 
problémes analogues sur les segments || Oz, avech,, = y(%o, Vo, 2). En intégrant 
et sommant toutes les inégalités du type (14) ci-dessus, nous obtenons, si la 


“4 
condition ener ee (15) 


est partout satisfaite, l’évaluation 


4a C = D(u) & [fod + 2am 4 uh + 2B w+ ut + 2yu) aV 
G 


(16) 
Bf] i+ 2at+ G+ 2Be+K + 29h av. 


9.1.3. Cette inégalité est d’un usage malaisé, car il faudrait encore résoudre 
tous les problémes auxiliaires 4 une dimension. De fagon analogue au § 3, 
retournons donc les problémes auxiliaires: partons de leur solution! Choisissons 
trois fonctions: f(x, y, z) continue en x et deux fois dérivable par rapport a »; 
g(x, y, z) continue en y, etc.; h(x, y, 2) continue en z, etc.; valant toutes ZérO 
sur J, et 1 sur I‘, et satisfaisant a la condition 


ie es guy + ee = 0 7 : (15’) 


-Alors nous choisissons 


01%, V2) = fags BU%) Ys 2) = Byys V2) = Mees 
et nous obtenons 


txt B ff B+ At lest B+ 28 ey the + 2h he) a, 
G 
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soit, en intégrant par parties, 


ZAMP 


of On Og Ov oh ‘ 02 — |] 2 2 2 dV 
4nC 224) (55° on te oy On =e Oz x) dS (Gee Bae 2 


sous la condition (15’). 


9.1.4. Définissons le champ vectoriel 


_ [or ee One 
P. Vox? Oy ae} 


la condition (15’) devient 
divp <0 


et notre évaluation (16’): 


G 


nC 22 Pin toe eae 


C’est le principe de THomson, appliqué ici aux champs p définis par (17) avec 


fo Per = Orsur 15) =". surl,: 


(16") 


(17) 


(15") 


(16") 


Pour obtenir de bonnes évaluations, on appliquera (16’) de préférence a 


des fonctions f, g, / satisfaisant méme f,, + g,, + 4,, = 0 dans tout G. 


9.2. Rigidité a la torsion 


9.2.1. G est un domaine plan simplement connexe, de contour J"; u(x, y) est 


P=Du) =2]] ud. 


déterminée par Au = — 2 dans G, w= 0 sur I’; la rigidit 


la torsion est 


9.2.2. Problémes auxiliaires 4 une dimension: choisissons deux fonctions 
a(x, y) et B(x, y); sur touc segment Tt || Ox dans G, a extrémités sur J’, résolvons 
le probléme aux limites: f(x, y)) = 0 aux extrémités, f,,(x, Yop) = a(%, Vo)3 f a 


une propriété extrémale, d’ot i (w2+ 2au)dx= iS 


(f2 + 2a f) dx. De méme, 


imposons g(%», y) = 0 sur J” et g,,(%9, y) = B(x, y). Intégrons et sommons les 


inégalités: 


[1+ max, (2+ A] PS | (+ 2af+e8+2B ead. 
G 


0.2.93 Retournons les problémes auxiliaires: choisissons f(x,y) continue en | 
x et deux fois dérivable par rapport a x; g(x, y) continue en y, etc.; toutes deux - 


s'annulant sur J” 
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Prenons alors «(x, y) = f,,, B(x, vy) = &,y; nous obtenons 


[b+ maxg (fee + Byy)] PS |] (E+ 2% het f+ 2g g,,) dd =— ff (P+ ead 
G G 
P S [ming(— fe — yy) — 11 ff (2 + @) a . (18) 


G 


Cette évaluation est optimalisée, parmi les multiples (kf, kg) du couple de 
fonctions (f, g) par le choix k = 2/ming (— /,, — £yy)3 dou 


[e+ 9) 24 
P<4— (18) 


ming = feos = Sux)? : 


- ou bien, en normant ming (— /,, — g,,) = 2: 


ap 


Si. —f.=#,, 22 dans G, alors P= // (7 + g2) dA. (18”) 


“ 


a 
On a deux fois l’égalité si f(x, y) = g(x, y) = u(x, ¥). 


9.2.4. Définissons le champ vectoriel 


p={—j5° — 2}; (19) 


nous avons: Si div p = 2 dans tout G, alors P < if p’ dA, c’est le principe de 
¢ 


THoMSON pour le probléme de la torsion [8], appliqué a nos champs p définis 
par (19) avec f= g=0Osur I. 
Pour obtenir de bonnes bornes supérieures, on appliquera (18") de préfé- 
rence a des fonctions f, g satisfaisant — /,, — g,, = 2 dans tout G. ; 
Remarque. En vertu de (18’), le principe de THomson peut aussi s’écrire 


2 — Max ming ie tow - Syy) 


choix de / et g Se 
J+ gaa 
G 


forme analogue a notre (7). 


: 
” 


~~ 


is ho | 


§ 10. Interprétation physique 


10.1. Considérons, pour simplifier, une membrane G fixée: u=0 sur le 
contour J". Le probléme de la fréquence fondamentale est généralement énoncé 
ainsi: 
«Quelle est la premiére valeur propre A, d'une membrane homogéne de masse 
-spécifique (par cm?) 1 sur le domaine G?» 


SH ee ee 
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Une question équivalente est la suivante: 

«Quelle est la masse spécifique A, d’une membrane homogéene sur G, vibrant a 
la fréquence fondamentale 1?» 

10.2. Interprétation des inégalités de Barta-Pélya. Soit (ct. § 3.3.3) (x, ¥) ae 
fonction positive continue et deux fois dérivable dans G, nulle sur T On saik 
alors que la membrane sur G, fixée sur son contour et de masse spécifique 
— Af|f, vibre a la fréquence fondamentale 1, la fonction propre correspondante 
étant (x, ¥). 

En effet, soit v(x, y) continue et dérivable dans G, et nulle sur J’, alors 

Af 


grad?v + 7 Ue = div( + grad ') + grad*v + a grad?/ — 2 “- grad v grad/ 


= div( + grad ') + (grade — = grad ') 


d’ou en intégrant: 


D(v) + fF vdA =0; 
H 


1s écifique = (— [vj = = D(v) 
se spécifique = (— A///)\” C4 ee 


tandis que Ryasse spécifique=(— at/pUf] = 1. 

Nous utilisons maintenant un principe général sur les systémes vibrants 
(CouRANT-HILBErRT I [1], p. 357): une diminution des masses ne peut 
qu’élever la fréquence, une augmentation des masses ne peut que l’abaisser ; par 
conséquent, une membrane homogene sur G, vibrant a la fréquence fondamentale 1, 
dot avo une densité de masse A, comprise entre 


IV 


1, 


Af Af 
it eae eee ee 
infg( 7 ) Shy S supe( 7 ie 
Ces inégalités ne sont autres que celles de BARTA-P6LyYa (cf. § 3. 3. S)e ‘ 
10.3. Nous considérons maintenant (cf. § 3.2) deux fonctions positives 
données: f(x, y) continue en x, etc.; g(x, y) continue en y, etc.; pour simplifier, 
supposons f = g = 0 sur I. : ‘ 
jeprsoupons un exemplaire du domaine G en laniéres horizontales de largeur 
dy; chacune sera assimilée A une corde vibrante de masse spécifique (par unité 
de longueur) = (—fealf) dy et de module d’élasticité (4 une dimension) = 1 dy; 
toutes ces cordes vibrantes ont la fréquence fondamentale 1 et la fonction 


propre correspondante f(x, yp») : 


feet (2) 720, 


de mé cee 
1éme, découpons un autre exemplaire de G en laniéres verticales de largeur 


Vol. XI, 1960 Sur la fréquence fondamentale d'une membrane vibrante 411 


dx; chacune sera assimilée A une corde vibrante de masse spécifique = (—g,,/g) dx 
et de module d’élasticité = 1 dx; ces cordes vibrantes, elles aussi, ont toutes 
la fréquence fondamentale 1. Nous avons donc un systéme de cordes indépen- 
dantes ; le systéme entier vibre avec la fréquence fondamentale 1. (La valeur propre 1 
a une dégénérescence -> co lorsque dx > 0 et dy > 0.) 


10.4. Modifions maintenant le systéme vibrant en lui imposant des con- 
traintes supplémentaires: nous «soudons» en chaque point (x, yo) la corde 
horizontale d’ordonnée yy et la corde verticale d’abscisse x9. Quelle est l’influence 
de ces contraintes supplémentaires sur la fréquence fondamentale? En vertu 
d’un autre principe général sur les systémes vibrants (cf. [1], p. 354), celle-ci 
ne peut en étre qu’augmentée. Nous obtenons donc une membrane (inhomogéne 
en général) de masse spécifique — f,,/f — g,,/g et de fréquence fondamentale = 1 
(le module d’élasticité étant égal a 1). 

Comme une diminution des masses ne peut provoquer qu’une nouvelle 
augmentation de la fréquence fondamentale ({1], p. 357), nous avons le 
théoréme: 

Toute membrane inhomogene de masse spécifique o(x, vy) S — f,,/f — g,,/g (et 
de module d’élasticité 1): 


Au+po(%,y)u=0 dans G e u=0 sw I, 


a une prenuere valeur propre u, = 1. 


10.5. Il en est ainsi pour la membrane homogéne de masse spécifique 
o = infg (— f../f — 8,,/g); donc: 
Une membrane homogéne sur G, vibrant a la fréquence fondamentale 1, dott 
avoir une masse spécifique A, = inf (— f,./f — &y,/g)- 
Cette interprétation est trés voisine du raisonnement qui sert de base a la 
«méthode de découpage» ([2], cf. aussi [4]). On n’oubliera cependant pas la 
-remarque 3.3.4: les cordes vibrantes auxiliaires peuvent avoir, par endroits, 
une masse spécifique négative. 


10.6. Notre principe se laisserait peut-étre généraliser par l’introduction 
d'une liaison élastique, en chaque point (%p, Yo), entre les deux cordes auxi- 
‘liaires qui s’y croisent: ce qui serait analogue aux «problémes intermédiaires » 
de WEINSTEIN et ARONSZAJN. De ce fait, la dégénérescence (cf. § 10.3) serait 
levée: le principe deviendrait applicable également aux valeurs propres supé- 
rieures. 


§ 11. Adaptation 4 l’équation de Schrodinger 
L’équation de SCHRGDINGER indépendante du temps est de la forme 


Au(x, y, 2) + [A— W(x, y, 2)] w(x, y, 2) = 0, 


PE ee NENT 
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ot W(x, y, 2) est une fonction donnée dans l’espace; en outre, on doit imposer 
a u(x, y, 2) des conditions aux limites, qui assurent l’intégration par parties © 


D(w) = [J] a — W) wav. 


Adaptons par exemple l’expression du principe donnée au § 3: Soient de 
nouveau /(x, y, z) > 0 continue en x et deux fois dérivable par rapport a Bp 
g(x, y, 2) > 0 continue en y, etc., h(x, y, 2) > O continue en z, etc. Sous véserve 
des conditions aux limites a imposer a f, g, h, on obtient a partir de (6): 


0 < D(u) + |] (@ at a al *e) uz dV 
= [f(r Ware lez an =e + “42 uw dV; 


G RE Suy hee 
A, = int[W(w, y 2) — 2. ae ae (20) 
Comme au § 7.1, on a l’égalité si f(x, y, z) = D(y, z) u(x, y, 2) 
g(x, y, z) = Pz, x) w(x, y, 2), W(x, y, 2) = Q(x, y) we, Y, 2). 
L’expression vectorielle du § 4 s’adapte ainsi, sous réserve des conditions aux 
limites: 


d’ou 


A, = inf[W(zx, y, 2) + div p — p*] . (20’) 


On a l’égalité si p = — grad u/u. ; 
On se reportera au § 4.3 pour les discontinuités admises du champ p. Une 
bonne borne inférieure sera obtenue si l’on définit p en sorte que 


W(x, y, 2) + div p — p? = constante. 
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Summary 


The starting point of the present paper is a simple method of PAYNE-WEINBER- 
GER to obtain lower bounds for the first eigenvalue of a membrane through com- 
parison with the Rayleigh principle for homogeneous strings. A generalization of 
this idea to auxiliary non-homogeneous strings leads to a maximum principle (7) 
(of which an inequality of BarTa-P6rva is a special case) or (7’); in its form (8), 
it was essentially known (although with some unnecessary restrictive hypotheses) 
to M. H. ProtTer in 1958. — This principle is closely related to the Thomson 
principle of boundary value problems; in particular, the allowed discontinuities of 
the concurrent vector fields (in the formulation (7’)) are the same as in THomsSoON’s 
principle. — Some applications of the principle to the calculation of rather sharp 
lower bounds are indicated, as well as an intuitive physical interpretation (in terms 
of variation of masses and constraints). It is valid in three dimensions as well; we 
indicate its extension to SCHRODINGER’S equation. 


Post-scriptum 


(a) Complément au § 4.3.3 et au $10: la formulation vectorielle du principe 
de maximum admet également une interprétation physique simple. 

(b) On aura avantage a lire aussi le travail, tres voisin, de M. H. PROTTER: 
Vibration of a nonhomogeneous membrane, Pacific J. Math. 9, 1249-1255 (1959). 


(Regu: le 6 février 1960.) 


Uber den Zusammenhang zwischen der Fliessbedingung 
eines starrplastischen Korpers und seinem Fliessgesetz 
Von Hans ZIEGLER, ETH, Ziirich 


1. Einleitung 


In der Plastizitatstheorie wird heute angenommen, dass zwischen der 


4 Fliessbedingung eines Kérpers und seinem Fliessgesetz eine enge Verkniipfung 


bestehe. Diese ist erstmals von R. v. Mises [1]) postuliert worden und wird 


_gewohnlich als Theorie des plastischen Potentials bezeichnet. Man kann sie, 


wenn man sich auf starrplastische Kérper beschrankt, in der von W. PRAGER 


_ ({2], S. 18) verallgemeinerten Gestalt, aber ohne Verwendung der Potential- 


funktion, etwa wie folgt formulieren: 


1) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 426. 


RATE 
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Es seien ¢, (k = 1,2, .--, n) die verallgemeinerten Verformungen des starr- 
plastischen Kérpers, das heisst ein Satz von Lagekoordinaten, welcher seine 
Deformation mit der gewiinschten Genauigkeit beschreibt. Ferner sei (unter 
Verwendung der bekannten Summationsregel) 


dA = Q, dq, (1.1) 


die elementare Deformations- bzw. Dissipationsarbeit. Dann kann man die 

QO, (k=1,2,...,n) als verallgemeinerte Krafte bzw. nach W. PRAGER als 

verallgemeinerte Spannungen bezeichnen. Sowohl der Verformungs- wie der 

Spannungszustand lassen sich durch Punkte mit den Fahrstrahlen gq bzw. Q— 
im n-dimensionalen euklidischen Spannungsraum R, darstellen. Definiert man 

hier zum Zwecke der Metrisierung das Skalarprodukt zweier Vektoren mit 


Q-q=Q%, | (1.2) 
so ist die elementare Dissipationsarbeit (1.1) durch 
dA =Q-dq (1.3). 


gegeben. Es entspricht dann jedem verallgemeinerten Verformungs- bzw. 
Spannungszustand q, bzw. Q, ein Vektor q bzw. Q oder, wie wir im folgenden — 
auch sagen werden, ein Punkt q bzw. Q in R, und umgekehrt. 

Betrachtet man den K6rper in einem gegebenen Stadium des Verformungs-_ 
prozesses, so lassen sich vom momentanen Spannungszustand aus gewisse 
andere ohne plastischen Fluss realisieren. Diese sollen als aplastische Spannungs- 
zustdnde und die Gesamtheit ihrer Bildpunkte in R, als aplastischer Bereich R,, 
bezeichnet werden. Ferner sei unter einem aplastischen Spannungsinkrement 
dQ* ein solches verstanden, das zwei benachbarte aplastische Spannungszu- 
stinde Q und Q’ ineinander iiberfiihrt. Die Fliessgrenze, eine Hyperflache in” 
R,, wird dann als Oberflache des aplastischen Bereiches durch diejenigen — 
aplastischen Punkte Q gebildet, fiir die nicht jedes hinreichend kleine Span-_ 
nungsinkrement dQ aplastisch ist. Liegen die benachbarten Punkte Q und Q’ 
an der Fliessgrenze, dann nennt man das aplastische Spannungsinkrement dQ* 
neutral. Liegt von den beiden Punkten genau einer an der Fliessgrenze, dann 
stellt dQ* eine Entlastung bzw. Belastung dar. 


Fiir die im aplastischen Bereich R,,, aber nicht an der Fliessgrenze liegenden_ 
Spannungspunkte ist das Verformungsinkrement dq null. Fiir die Punkte an_ 


der Fliessgrenze kann es gleich oder ungleich null sein. Beziiglich seiner Rich- 
tung gelten folgende Sitze: 


: 
Satz 1: Mit Q’ und Q” ist auch jeder Spannungszustand } 
O= 0 +4.(0" =O), (Oneaaem al 
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aplastisch. Mit anderen Worten: Mit zwei beliebigen Punkten Q’ und Q" liegt 
auch ihre ganze Verbindungsstrecke im aplastischen Bereich, das heisst der 
aplastische Bereich ist konvex. 


Satz 2: Ist Q ein Spannungszustand an der Fliessgrenze und dQ* ein be- 
lhebiges aplastisches Spannungsinkrement, das vom wirklichen Spannungsin- 
krement dQ véllig unabhangig ist, so gilt 


dQ* -dq <0. (1.5) 


Mit andern Worten: Die Normalprojektion des Verformungsinkrements auf 
jedes mit dem wirklichen Spannungszustand vertragliche, aplastische Span- 
nungsinkrement ist nichtpositiv. 


Satz 3: Fiir das wirkliche Spannungsinkrement dQ gilt 
dQ -dq=0. (1.6) 


Das heisst: Die Normalprojektion des Verzerrungsinkrements auf das wirk- 
liche Spannungsinkrement ist nichtnegativ. . 

Die hier gewahlte Gestalt der Satze 1 bis 3 ist mit einer Formulierung von 
D. C. DRUCKER [3] aquivalent. Sie vermeidet die Beniitzung der Fliessflache 
mit ihren allfalligen Singularitaten und ins Unendliche reichenden Teilen und 


lasst sich zudem leicht auf den Funktionenraum iibertragen. Wird R, in Q 
durch ein regulares Stiick der Fliessflache begrenzt, so folgt aus (1.5), dass dq in 
ihre aussere Normale fallt. Hat die Fliessflache in Q eine Kante oder Ecke, so 
kann diese nach (1.4) nur ausspringend sein; dq wird durch (1.5) auf den 
Winkelraum beschrénkt, der durch die Grenzlagen der ausseren Normale bei 
der Annaherung an Q gebildet wird, und schliesst zudem nach (1.6) mit dQ 
einen nichtstumpfen Winkel ein. 


2. Problemstellung 


R. v. Mises hat den durch die Satze 1 bis 3 beschriebenen Sachverhalt 
anders ausgesprochen. Hervorzuheben ist vor allem, dass er ihn ausdriicklich 
_ als Hypothese und zudem unter Beschrankung auf ein Raumelement statuiert 
_ hat, indem er als Verformungen g, und Spannurgen Q, die lokalen Verzerrungs- 

komponenten ¢;; und die Spannungskomponenten o,, benitzte. Es ist das 
- Verdienst von W. PRAGER, die Theorie im Sinne von Abschnitt 1 verallgemei- 
-nert zu haben. Gegeniiber der y. Misesschen Fassung bedeutet dies eine 
- wesentliche Erweiterung, deren praktische Bedeutung auf der Hand liegt; am 
3 hypothetischen Charakter der Theorie andert aber diese Neufassung nichts. 
| Es sind bisher drei Versuche unternommen worden, die Satze 1 bis 3 zu 
- begriinden. Zundchst haben J. F. W. Brsuor und R. H1t [4] gezeigt, dass die 
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Satze 1 und 2 fiir das polykristalline Raumelement aus einer plausiblen Hypo- 
these iiber das Verhalten des Einkristalls beim plastischen Gleiten folgen. 
Sodann hat D. C. DRUCKER [3] nachgewiesen, dass sie in der v. Misesschen 
Fassung richtig sind, wenn man iiber das Verhalten des Raumelements zwei 
Voraussetzungen machen darf. Er postuliert, dass dann, wenn man am be- 
liebig vorbelasteten Element eine zusatzliche Belastung langsam anbringt und 
wieder entfernt, die durch diese Zusatzlasten geleistete Arbeit wahrend der 
Belastung positiv und tiber den ganzen Prozess nichtnegativ sel. Diese Voraus- 
setzungen sind indessen nicht ganz so selbstverstandlich, wie sie scheinen 
mégen. So ist zum Beispiel ein Element aus einem Material mit dem Span- 
nungs-Dehnungs-Diagramm von Figur 1, wenn es in einachsigem Zug bis an 


o 


% @ 
A B 


E 
Figur 1 
Spannungs-Dehnungs-Diagramm eines starrplastischen Materials mit Verfestigung. 


die Fliessgrenze o, vorbelastet wird, unter konstant gehaltener Belastung 
im labilen Gleichgewicht, da sich mit dem Beginn des Fliessens infolge der — 
Querkontraktion die Normalspannung erhéht. Um den Vorgang quasistationar 
zu halten, muss eine zusatzliche Druckkraft eingefiithrt werden, und wenn die 
bei B einsetzende Verfestigung geniigend stark ist, gibt es auf dem zugeh6rigen — 
Teil des Diagramms einen Punkt C, in dem die Zusatzlast wieder auf null ab-— 
sinkt. Ihre Arbeit ist sowohl wahrend des Aufbringens wie tiber den ganzen 
Zyklus negativ und erfiillt damit keines der beiden Postulate. f 
Ein Versuch ganz anderer Art, die Satze 1 und 2 zu begriinden, wurde vom 
Verfasser [5] unternommen, und zwar im Rahmen einer Erweiterung der 
Theorie irreversibler Prozesse von L. ONSAGER [6]. Auch diese Beweisfiihrung © 
beruht auf zwei Postulaten, namlich auf der Annahme, dass sich irreversible 
Vorgange unter mikroskopischen Schwankungen abspielen, die erstens selbst 
reversibel und zweitens im Vergleich zu den Molekularbewegungen langsam — 
sind. Das sind im wesentlichen die Voraussetzungen, die schon der Onsager- 
schen Theorie zugrunde liegen. Aus ihnen kann meines Erachtens auf die 
Richtigkeit der Satze 1 und 2 direkt in ihrer verallgemeinerten Pragerschen | 
Form geschlossen werden. Die in [6] angewandte Schlussweise macht aber von 
statistisch-thermodynamischen Anschauungen Gebrauch und liegt damit aus- 
serhalb der Mechanik im engeren Sinn. Sie ist zudem zweifellos im einzelnen 
noch verbesserungsfahig. Tatsache ist jedenfalls, dass sie zurzeit noch nicht 
als anerkannt gelten kann, auch wenn sie kaum auf Kritik gestossen ist. 
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In dieser Situation scheint es nicht unangebracht, sich fiir einmal zu be- 
scheiden und unter Verzicht auf alle physikalischen Postulate wenigstens in 
aller Strenge nachzuweisen, dass aus “der Annahme, die Sitze 1 bis 3 seien in 
der v. Misesschen Fassung, das heisst fiir das Raumelement richtig, auch ihre 
Giiltigkeit in der Pragerschen Verallgemeinerung, das heisst im Sinne von 
Abschnitt 1 folgt. Dieser Nachweis soll im folgenden gebracht werden, und 
zwar unter Verwendung einer neulich von J. L. SyNcr [7] in anderem Zusam- 
menhang beschriebenen Darstellungsart. 


3. Ein Beispiel 


Zunichst soll ein Beispiel, das von W. PRAGER ([2], S. 44) diskutiert 
worden ist, die Satze 1 bis 3 illustrieren und gleichzeitig als Testfall fiir 
die weitere Orientierung dienen. Fiir ein einzelnes Raumelement kann man 
als Verformungen die Hauptdehnungen ¢, ¢,,¢, und als Spannungen die 
Hauptspannungen oj, 02,0, verwenden, da dA = a; de; die elementare Dissi- 
pationsarbeit ist. Unter der bekannten Fliessbedingung von R. v. MIsEs ist der 


aplastische Bereich R, ein Kreiszylinder vom Radius (2/3)1/2o, wenn oy die 


_ Fhessgrenze fiir den linearen Spannungszustand bezeichnet, und die Achse 


dieses Zylinders bildet mit den positiven Koordinatenachsen in R, gleiche 
Winkel. Unter der Fliessbedingung von H. Tresca dagegen ist R, ein regulares 
Prisma mit hexagonalem Querschnitt, das dem v. Misesschen Zylinder einge- 
schrieben ist. Figur 2 zeigt beide K6rper, in axialer Richtung betrachtet, und 


_ Figur 2 
Aplastische Bereiche nach R. v. Mises bzw. H. Tresca im Hauptspannungsraum. 


zwar unter der Annahme, dass noch keine Verfestigung stattgefunden habe. 
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Satz 1 ist offensichtlich in beiden Fallen erfiillt, und nach Satz 2 ist der Vektor 
de, déy, de, des Verzerrungsinkrements entweder normal zur Oberflache nach 
aussen gerichtet, oder er ist normal zu einer Kante, liegt 1m ausseren Winkel- 
raum, der durch die Normalebenen durch diese Kante zu den anschliessenden 


Seiten gebildet wird, ‘und zwar so, dass der Winkel zwischen ihm und dem 
Vektor do,, doz, do, nicht stumpf ist. 
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Es sei nun ein Balken mit rechteckigem Querschnitt aus diesem von 
jetzt an als idealplastisch vorausgesetzten Material auf Zug und Biegung 
beansprucht und hievon in Figur 3 eine Scheibe der Dicke 1 betrachtet. 


Figur 3 
Elementarscheibe des auf Zug und Biegung beanspruchten Balkens. 


Der Verformungs- und der Spannungszustand in dieser Scheibe sind in 
Wirklichkeit kompliziert und werden, falls man die Langeneinheit hinrei- 
chend klein gewahlt hat, durch je sechs von y und z abhangige Verzerrungs- 


bzw. Spannungskomponenten beschrieben. Man kann aber das Problem im ~ 


Sinne der Festigkeitslehre durch die Annahme vereinfachen, dass die Quer- 


schnitte bei der Verformung eben bleiben und sich wie starre Kérper bewegen. — 
Um mit der Inkompressibilitat des Materials nicht in Konflikt zu kommen, | 


muss man sich dabei freilich vorstellen, dass der Zusammenhang zwischen den 
Langsfasern gelést sei, und das ist natiirlich nur dann erlaubt, wenn der Span- 


nungszustand im wesentlichen linear mit x als Achse ist. Wie Figur 4 zeigt, — 


wird die Verformung der Scheibe in dieser Naherung durch die Dehnung ¢ sowie 
die Kriimmung der elastischen Linie, namlich den Winkel x zwischen den 


beiden Stirnflachen, beschrieben, und die zugehdrigen verallgemeinerten Span- — 


nungen sind die Normalkraft N sowie das Biegemoment M, da die elementare 
Dissipationsarbeit durch dA = N de + M dx gegeben ist. 


f 
Iz 
qe 
Figur 4 
Verformung der Scheibe. 


Figur 5 
Plastische Spannungsverteilung in der Scheibe. 


Da das Material als starrplastisch vorausgesetzt ist, kommt es zur Verfor- 
mung gemiss Figur 4 erst dann, wenn der ganze Querschnitt plastisch geworden 
ist. Die Lage der neutralen Faserschicht ist durch 1 gegeben, und die Normal- 
spannung ist bei Beginn des Flusses gemiss Figur 5 iiber den Querschnitt 


( 


cee ae eae 


s 
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verteilt. Dabei gilt, wie man leicht tiberpriift, 


Et x 
und (3.1) 
Line i aye cal gate “ h _ ] 
N=+2bro,, Ms hee 4 #2) 9, (-3 srs (S22) 
bzw. 
i ta Me ah h 

t+bho,, M=0 (rs eo r2>)- (353) 

Der aplastische Bereich R, wird also (Figur 6) durch die beiden Parabein 

bh? N \2 
M=-+- — |. .— 

Lp 01 (aera) | i222) 


begrenzt. Er ist konvex, und man tiberzeugt sich leicht davon, dass das Ver- 
formungsinkrement de, dx in die aussere Normale bzw. in einen der beiden 
schraffierten Winkelraume, beim Abwandern des Spannungspunktes aus einer 
Ecke aber auf den naheren Rand des zugehérigen Winkelraumes fallt, so dass 
auch in dieser Darstellung die Satze 1 bis 3 gelten. 


M 


Figur 6 
Aplastischer Bereich fiir die auf Zug und Biegung beanspruchte Scheibe. 


Wie in diesem Beispiel, so wird auch in jedem anderen Fall der Ubergang 
vom wirklichen Problem zur vereinfachten Darstellung gemass Abschnitt 1 
dadurch vollzogen, dass man den Freiheitsgrad des K6rpers beschrankt und 
nur Verformungen in Betracht zieht, die sich durch wenige Parameter 9, 


beschreiben lassen. Diese speziellen Verformungen brauchen nicht kinematisch 
_ zulassig, das heisst mit den kinematischen Randbedingungen, den Kontinui- 
tatsbedingungen und allenfalls der Forderung der Inkompressibilitat vertrag- 
lich zu sein, auch wenn die Vereinfachung natiirlich um so brauchbarer ist, je 


_ 


_ weniger von diesen Bedingungen verletzt werden und je besser die zugehérigen 
- Spannungszustinde mit den wirklichen iibereinstimmen. Die verallgemeinerten 


_ Spannungen Q,, welche den Verformungen 4;, entsprechen, ergeben sich aus der 
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Darstellung (1.1) der elementaren Dissipationsarbeit und werden im allge- 
meinen nicht statisch zulassig, das heisst mit allen dynamischen Randbedin- 
gungen sowie den Gleichgewichtsforderungen im Einklang sein. 

Um den aplastischen Bereich R, zu erhalten, muss man diejenigen Span- 
nungszustande im Kérper angeben konnen, welche fiir verschiedene Verhalt- 
nisse der Q, erstens an der Fliessgrenze liegen und zweitens mit den gewahlten 
speziellen Verformungszustaénden vertraglich sind. Das ist um so leichter, je 
einfacher die durch die g, dargestellten Verformungszustande sind; auf alle 
Falle steht und fallt mit dieser Méglichkeit die Brauchbarkeit der Vereinfa- 
chung. Es mag gelegentlich von Vorteil sein, wenn man den K6rper als ela- 
stisch-plastisch betrachtet und von der elastischen Spannungsverteilung aus 
auf diejenige an der Fliessgrenze schliesst. Fiir das Folgende ist indessen nur 
die Spannungsverteilung an der Fliessgrenze von Bedeutung, und es andert 
sich an den Ergebnissen nichts, wenn man sich vorstellt, dass bei gegebenem 
Verhaltnis zwischen den Q, der Spannungsverlauf unterhalb und an der Fliess- 
grenze nur in der Intensitat, nicht aber in der Form verschieden sei. 


Der plastische Fluss setzt im besprochenen Beispiel erst dann ein, wenn ~ 


alle Elemente der Scheibe an der Fliessgrenze angelangt sind. Im allgemeinen 


ist dies nicht nétig. Wenn etwa der ganze Balken betrachtet und angenommen ~ 
wird, dass nur die Scheibe von Figur 3 gefahrdet sei, dann liegt ein Beispiel — 


fiir den Fall vor, wo der K6rper als solcher die Fliessgrenze erreicht, ohne dass 
dies alle seine Elemente tun. 


4. Der lokale Spannungszustand 


Um den am Schluss von Abschnitt 2 in Aussicht gestellten Beweis zu — 


erbringen, seien die Verzerrungs- und Spannungskomponenten in einem belie- 
bigen Element des betrachteten Kérpers mit e,; bzw. o;; bezeichnet und durch 
die Vektoren e bzw. s im neundimensionalen euklidischen Spannungsraum Ry, 
dargestellt. Definiert man hier das Skalarprodukt mit 


S°e= Oi; Ei» (4.1) 
so ist die elementare Dissipationsarbeit je Raumeinheit durch 


gegeben pr — 


Mit Riicksicht auf die Symmetrie der Tensoren e;; und o,, liegen die Vekto- 


ren e und s in Wirklichkeit in einem sechsdimensionalen linearen Unterraum — 


s von R,. Hier entspricht jedem Verzerrungs- bzw. Spannungszustand e, j 
ZW. 0; 


LW. Oi; ein Vektor e baw. s und umgekehrt, und der aplastische Bereich wird 
wie in Abschnitt 1 definiert. Er kann sich wie in den Fallen von Figur 2 ins 


Unendliche erstrecken , So schliesst man ja gewodhnlich plastischen Fluss unter 
einem hydrostatischen Spannungszustand aus, 


: 
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Es sei von jetzt an vorausgesetzt, dass fiir den lokalen Spannungszustand 
und das zugehérige Verzerrungsinkrement die Satze 1 bis 3 gelten. Dann ist 
erstens nach (1.4) mit s’ und s” auch jeder Spannungszustand 


s=s'+a(s"—s'), (0<a<1) 
oder in anderer Schreibweise mit o;; und oj; auch jeder Spannungszustand 
G,;=9;,+a(ci,-—9;), (0<a<1) (4.3) 


aplastisch. Zweitens gilt nach (1.5), wenn s bzw. o;,; ein Spannungszustand an 
der Fliessgrenze und ds* bzw. do*, ein beliebiges aplastisches Spannungsinkre- 
ment ist, fiir das wirkliche Verzerrungsinkrement de bzw. de, ; die Ungleichung 
ds* - de < 0 oder 

dox de;, SO. (4.4) 


Drittens gilt nach (1.6) fiir das wirkliche Spannungsinkrement ds bzw. de die 
Ungleichung ds - de = 0 bzw. : 
do;,de,;,; 20. (4.5) 


Diese Voraussetzungen sind insbesondere in den beiden Fallen erfiillt, die 
in Abschnitt 3 im Zusammenhang mit Figur 2 diskutiert worden sind. 


5. Der globale Spannungszustand 


Der Verformungs- und der Spannungszustand des ganzen Kérpers werden 
durch die von den Koordinaten x, y und z abhangigen Funktionen é;; bzw. o;; 
dargestellt, und diese lassen sich [7] als Vektoren E bzw. S im Funktionen- 
raum F interpretieren. Da die elementare Dissipationsarbeit durch das tiber 
den ganzen Kérper K erstreckte Raumintegral 


As fous de, , dv (5.1) 
K 


gegeben ist, empfiehlt es sich, fir die Metrisierung des Funktionenraums das 
_ Skalarprodukt hier durch 
See | Oj 4 a (5.2) 


K 


zu definieren. Diese Definition ist zulassig (vgl. [7], S. 37), da ihr zufolge, wie 
man leicht iiberpriift, das kommutative Gesetz S:E=E-S, das distributive 
 S-(E+E')=S-E+S-E’, die Beziehung S-O =0 fir die Multiplikation mit 
~ dem Nullvektor und schliesslich das assoziative Gesetz (4 S)-E=a(S-E) fiir 
die Multiplikation mit einem Skalar gilt. Die Metrik in F ist dann positiv 


Li i 
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definit, und die elementare Dissipationsarbeit (5.1) ist durch 


dA=S:dE (5.3) 9 
gegeben. 

Jedem Verformungs- bzw. Spannungszustand im Kérper K entspricht ein 
Vektor im Funktionenraum. Dabei beschrankt sich diese Darstellung keines- 
wegs auf kinematisch zulassige Verformungs- und statisch zulassige Spannungs- 
zustande: vielmehr sollen alle denkbaren Zusténde zugelassen sein, fiir welche 


die Integrale (5.1) und (5.2) existieren. Der aplastische Bereich # im Funktio- 
nenraum wird analog wie in Abschnitt 1 definiert. Er kann sich ins Unendliche 
erstrecken; so nimmt man ja gewohnlich an, dass kein plastischer Fluss auf- 
trete, wenn der Spannungszustand iiberall in K hydrostatisch ist. Im allge- 
meinen werden gewisse Elemente des Kérpers K schon fiir Spannungszustande 
S plastisch, die noch nicht an der Fliessgrenze liegen. Plastischer Fluss ist aber 
erst médglich, wenn ein geniigend grosses Teilgebiet von K plastifiziert ist, und 
bei Fliessbeginn befinden sich im allgemeinen die Elemente von K in gewissen 
Teilbereichen an ihrer Fliessgrenze, in anderen noch darunter. 

Es sei nun angenommen, dass Satz 1 zwar fiir jedes Element von K, aber 
nicht fiir den ganzen Spannungszustand gelte. Dann gibt es zwei aplastische 
Spannungszustinde S$’ und S” sowie einen Wert a so, dass der Spannungs- — 
zustand 

S=S’'+a(S"—S'), 0<a<l) (5.4) @ 


nicht aplastisch ist. Das ist aber nur méglich, wenn in wenigstens einem Ele- 
ment von K die Fliessgrenze tiberschritten, mithin im Widerspruch zur ersten 
Annahme von Abschnitt 4 der Spannungszustand (4.3) nicht aplastisch ist, 
obschon es die Spannungszustande o;,; und oj; sind. 

Ebenso einfach ist der Nachweis, dass fiir einen Spannungszustand S$ an 
der Fliessgrenze, ein beliebiges aplastisches Spannungsinkrement d$* und das — 
wirkliche Verformungsinkrement dE die Ungleichung | 


dS* -dE <0 (5.5) 8 
bzw. nach (5.2) 


| dot de,, dv <0 (5.6) 
Ie 


besteht, falls fiir jedes Element von K Satz 2 richtig ist. Man kann namlich ~ 
diese Elemente in zwei Gruppen einteilen, je nachdem sie unter dem wirklichen 
Spannungsinkrement dS fliessen oder nicht. Fiir die Elemente der zweiten 
Gruppe ist de,; = 0; sie liefern also keine Beitrage zum Integral (5.6). Fiir . 
diejenigen der ersten Gruppe gilt (4.4), da dS* nur mit allen do* aplastisch ist ; / 
ihre Beitrage zu (5.6) sind demnach nichtpositiv. c 


| 
| 
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Vol. XI, 1960 Zusammenhang zwischen der Fliessbedingung eines Kérpers 423 


Schliesslich folgt aus der Giiltigkeit von Satz 3 fiir jedes Element von K 
dass fiir das wirkliche Spannungsinkrement d§ die Ungleichung 


aS -dE =O (5..4) 
bzw. 


[4o:; de,, dv =0 (5.8) 


k 


besteht. Die unter dS nichtfliessenden Elemente leisten namlich keine, die 
librigen nach (4.5) nichtnegative Beitrage zum Integral (5.8). 
Somit hat man 


Satz 4: Gelten die Satze 1 bis 3 fiir das Raumelement, dann sind sie auch 
fiir die Darstellung des ganzen Verformungs- und Spannungszustandes im 
Funktionenraum richtig. 


6. Behandlung im linearen Unterraum 


Wie in Abschnitt 3 ausgefiihrt wurde, bedeutet der Ubergang vom wirklichen 
Problem zur vereinfachten Darstellung von Abschnitt 1 die Beschrankung auf 
Verzerrungszustande, die sich durch Uberlagerung aus ” speziellen, voneinander 
unabhangigen Zustanden gewinnen lassen. Im Funktionenraum bedeutet dies 
stets die Restriktion auf einen linearen Unterraum F, (vgl. [7], S. 28), der den 
Ursprung enthalt und durch eine orthonormale Basis I, (k = 1, 2,... , m) dar- 
gestellt werden kann, so dass 

R= (E-I,) I, g(G.ck.) 


ist. Die zugehérigen Spannungszustande werden im allgemeinen nicht durch 
Vektoren § aus F, dargestellt. Indessen kann ein beliebiger Spannungszustand 
S (vgl. [7], S. 61) gemass 


Ge Ses SS 1 (6.2) 


in seine Normalprojektion S’ auf F, und einen zu F, orthogonalen Vektor Sp 
zerlegt werden, und da nach (5.3), (6.1) und (6.2) die elementare Dissipations- 


~ arbeit durch 


dA =S-dE=(S-I,)(GE-1,)=S'-dE (6.3) 


gegeben ist, liegt die Vermutung nahe, dass dennoch die Behandlung im linearen 


-Unterraum F, méglich sei. 
Diejenigen Vektoren S$’, die durch Normalprojektion aus aplastischen 


4 Spannungszustanden S hervorgehen, bilden die Normalprojektion F’ des apla- 
4 stischen Bereiches F auf F,. Ein Punkt, der F’ nicht angehért, kann nicht durch 


_ Projektion eines Punktes aus F erhalten werden. Somit gilt 
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Satz 5: Die Normalprojektion dS*’ jedes aplastischen Spannungsinkre- 
ments dS* liegt im Bereich P’. 

Ferner ist nach (6.2) die Normalprojektion des Vektors S = S,+a (S,—S)) 
durch S$’ = Si + a (S5 — S{) gegeben. Daher gilt 


Satz 6: Mit dem aplastischen Bereich F ist auch seine Normalprojektion F' 


auf den linearen Unterraum F, konvex. 
Die Beschrankung auf die speziellen Verformungszustande (6.1) hat zur 


Folge, dass fiir plastischen Fluss nur diejenigen Punkte S auf der Fliessgrenze — 


in F, das heisst auf der Oberflache von F, in Betracht zu ziehen sind, fiir die 
mindestens ein Verformungsinkrement dE dem linearen Unterraum F, ange- — 


hért, mithin von der Gestalt 


dE = (dE: I,) I; (6.4) — 


ist. Stellt dann dS* ein beliebiges aplastisches und dS das wirkliche Spannungs- 
inkrement dar, so gilt analog zu (6.3) 


dS*.-dE=dS*’-dE, dS -dE=dS'-dE (6.5) 


und nach (5.5) bzw. (5.7) 


dS*'.-dE<0. dS’-dE=0. (6.6) 


Aus der ersten Ungleichung (6.6) folgt zunachst, wenn man noch Satz 5 


beachtet: 


Satz 7: Diejenigen Spannungszustande, welche an der Fliessgrenze liegen — 
und mit Verformungszustanden der Gestalt (6.1) vertraglich sind, ergeben bei — 


der Normalprojektion auf F, die Randpunkte des Bereiches F’. 
Damit ist die Méglichkeit der Behandlung im Unterraum F, nachgewiesen, 


und die Satze 6 und 7 sowie die Ungleichungen (6.6) kénnen schliesslich wie — 


folgt zusammengefasst werden: 


Satz 8: Bei Beschrankung der Verformungszustande E im linearen Unter- — 


raum /,, von F kann man ganz in diesem Unterraum arbeiten, wenn man die 


Spannungszustande § und den aplastischen Bereich F durch ihre Normalpro- 


jektionen S’ baw. F’ ersetzt. In den projizierten Gréssen gelten die Satze 1 bis 3. 


7. Der verallgemeinerte Spannungszustand 


Wir verfiigen jetzt hinsichtlich der verallgemeinerten Verformungszustande 
von Abschnitt 1 tiber zwei Darstellungen, namlich durch Punkte E in F 
sowie durch Punkte q in R,. Ferner gehért zu jedem Spannungszustand, wel- 
cher mit dieser Einschrankung der Verformung vertraglich ist, ein Bildpunkt $’ 


tt a 
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in /, und ein solcher Q in R,. Diese Bildpunkte stellen nicht die ganzen 
Spannungszustande dar, nach (1.3) und (6.3) aber diejenigen Anteile hievon, 
die bei Verformungen der hier betrachteten speziellen Gestalt Arbeit zu leisten 
vermogen. Im Hinblick auf die Eindeutigkeit der virtuellen Arbeit entspricht 
einem gegebenen Spannungszustand nur je ein einziger Bildpunkt in den 
Spannungsraumen f, und R,. Zwischen den Punkten S$’ und Q besteht also 
eine umkehrbar eindeutige Zuordnung. Da die beiden Punkte den gleichen 
Spannungszustand darstellen, miissen zudem beide gleichzeitig innerhalb, aus- 
serhalb oder auf der Oberfliache des zugehérigen aplastischen Bereiches F’ 
baw. R, liegen. Damit entspricht insbesondere einem aplastischen Spannungs- 
inkrement dS*’ in F, ein solches dQ* in R,,. 

Da das Skalarprodukt in beiden Raumen F, und R, durch die Deformations- 
arbeit definiert ist, hat es fiir entsprechende Vektorpaare $’, dE bzw. Q, dq den- 
selben Wert. Somit folgt aus (6.6) die Richtigkeit von Satz 2 und 3 auch in der 
urspriinglichen Form (1.5) bzw. (1.6) von Abschnitt 1. Nimmt man ferner an, dass 
in F, mit S’ und $” auch jeder Punkt (5.4) im Inneren der Verbindungsstrecke 
von S$’ nach §” aplastisch sei, so folgt dies sofort auch fiir die Bildpunkte Q’, Q” 
und Q in R,. Um schliesslich noch zu zeigen, dass Q im Inneren der Strecke 
von Q’ nach Q” liegt, beachten wir, dass man die Ortsbeziehung (5.4) zwischen 
S’, S” und S bei geeigneter Wahl der orthonormalen Basis I, in F, auch durch 
die Feststellung ausdriicken kann, dass 


(S’—S)-,>0, (S—S')-E>0 (7.1) 


sei, wahrend die entsprechenden, mit den tibrigen Einheitsvektoren gebildeten 
Skalarprodukte samtlich verschwinden. Mit Riicksicht auf die Ubereinstim- 


- mung des Skalarproduktes in den Réumen F, und R, geht die orthonormale 


Basis I, von F, in eine solche J, von R,, iiber, und es folgt aus (7.1), dass 
(OeetO( I, 0, (FOS, 0 (7.2) 


ist, wahrend die entsprechenden Produkte mit den iibrigen Einheitsvektoren 
verschwinden. Damit ist aber bewiesen, dass mit F’ auch R, konvex ist, und 


Bes kommt daher schliesslich 


Satz 9: Gelten die Satze 1 und 2 fiir das Raumelement, dann sind sie auch 


- fiir die Darstellung der verallgemeinerten Verformungen und Spannungen in 


<4 
mi 


: 


Lan ee 


R,, richtig. 


Abstract 


This paper supplies a proof, free of any physical assumptions, that v. MisEs’ 
theory of the plastic potential, if valid tor the elements of a rigid-plastic body, also 


holds for the representation in generalized strains and stresses. 
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Laminar-Film Condensation on a Flat Plate in the Absence 
of a Body Force : 


By Rosert D. Cess, Pittsburgh 35, Pa., U.S.A.1) 


1. Introduction 


Analyses of laminar-film condensation generally deal with condensation of a * 
saturated vapor on a surface located within a gravitational field and with negligible — 
vapor velocity. This problem was originally studied by NussEett [1]*) for the case 
of a vertical flat plate, and both fluid acceleration and thermal convection within ~ 
the liquid film were neglected. SPARROW and GREGG [2] have recently obtained a_ 
solution to this problem without neglecting acceleration and convection effects. — 
This was accomplished through application of boundary-layer techniques and : 
subsequent utilization of a similarity transformation. . 

A further consideration of laminar-film condensation on vertical surfaces in- 
volves a finite vapor velocity. Here the motion of the liquid film is affected by both i 
the gravitational body force and the sweeping effect of the vapor. In this situation, 
a similarity transformation does not exist, and a complete solution to the problem 
would b2 quite complicated. Only approximate solutions have been obtained based 


upon the assumption that the shear stress at the liquid-vapor interface is a constant 
known quantity’). 


e 1) Heat Transfer and Flow Section, Mechanics Department, Westinghouse Research Labora- 
ories. 

*) Numbers in brackets refer to References, page 433. 

*) For example, see Reference [3] 
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The purpose of the present paper is to investigate laminar-film condensation 
on a flat plate for the case in which no body force is present. The motion of the 
liquid film therefore results solely from. the sweeping effect of the adjacent vapor. 
This problem is illustrated in Figure 1. A saturated vapor at temperature 7. flows 


Figure 1 


Physical model and coordinate system. 


with a free-stream velocity U,, across a cooled flat plate having a constant surface 
temperature 7,,. The vapor condenses on the plate surface producing a liquid film 
of thickness 0. 

Physically, such a situation can occur in one of two ways. First, if the plate is 
oriented normal to any acting body force; for example, if the plate is placed 
horizontally in a gravitational field. The second situation involves the complete 
absence of a body force as would occur in a nonrotating space vehicle. 


2. Basic Equations 


It will be assumed that a continuous laminar film of condensate exists. Let u 
and uv be the velocity components of the liquid, and U and V be those of the vapor. 
Temperature is denoted by T. For y < 6 the equations expressing the conservation 
of mass, momentum, and energy may be written as 


Ou Ov 

ae Sere a 1 
Tote Qi (1) 

Ou Ou 02u 
aS 2 
uU a Sat dy 1 oy?” (2) 

oT oT oT 
prs pet eee 3} 
nary ‘teal Ov " OV (3) 


while for the vapor (y > 6) the continuity and‘momentum equations become 


oU OV 
——— =) 4 
= >. ; (4) 
oU 0U 02U 
cated ae east 5 
U ee -b ay Vy dy? (5) 


_ The quantities « and v are the thermal diffusivity and kinematic viscosity, respec- 
_ tively. Further, the nomenclature will be adopted that a physical property with no 
subscript refers to the liquid, whereas the subscript v denotes the vapor. 
The boundary conditions at the plate surface are wu = v = 0 ands Dl == Ie while 
T = Tx. for y = 6 and U = U, for y = oo. Three additional boundary conditions 
are necessary. Across the liquid-vapor interface both velocity and shear stress 
must be continuous, so that within the iramework of boundary-layer theory u = U 
and pu du/dy = uw, OU/dy for y = 6, where mu is the dynamic viscosity. In addition, 


i rc 2 Oe fe 


Acs | 
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continuity of mass across y = 6 requires that 


do dé 
hige eat ee v) 
: (« dx hey ey ( dx y=6 


with 9 denoting density. 
Equations (1) and (4) are satisfied by introducing stream functions py and W 
defined by the relations 


oy __ oy eae oa 
iC imaeoara Uo tye 5) = ov , Ox ~ 
Defining further 
ys LL, (oo) 
a ae F(t = O(n) = =, nav , 
f(n) V0, Uno x , (7) %» U x ee — ess, VV, x 


equations (2), (5), and (3) transform, respectively, to the ordinary differential — 
equations 


mt 1 Vy We 6 
f+ (2) ti =0, (6) 
1 
Fu a a jE ie = F (7) | 
Q" 4 So (lee =, (8) 
2 v } 
The boundary conditions become 

f(0) = #(0) = 0, (0) =1, (9a) 


fins) = (42) Few) » 11a) = Fema), 1") = (22%) FM), Ole) = 9, (90) 


~~ e (co) ae (9c) 
where 75 = 6 //U/v, % - 

One may note that equations (7) and (9b) are analogous to vectored blowing or — 
suction, where in this case the vectored blowing or suction is specified at 7 = n5_ 
rather than 7 = 0. In view of this result, consider a function G(n) defined by ; 

if , 
GU + a GG—70 (10) 
with the boundary conditions 


G(0) =G,,, G’(0) = 0, Goo) = 1. ; 


This corresponds to the problem of normal blowing or suction as treated by Em- t 
MONS and LEIGH [4], and it has been shown that vectored blowing or suction can _ 
be determined from the solution of equation (10) with its boundary conditions — 
through use of a fictitious origin for 7. Consequently, the solution of equation (7) 
may be written as } 
F(n) =G (q — »), (11) 


f 
where y and the appropriate value for G,, are constants to be determined from nel 
boundary conditions at 7 = np. ; 


- 
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The quantity 1,» remains to be evaluated, and this may be accomplished in the 
same manner as considered by SPARROW and GreEaG [2]. An over-all energy balance 
on the liquid film yields 

é 0) 


p 


k : oT - ‘ TT 
al Se): A dx#=0 hry | udy + @ e» | u (To — T) dy, 
0 


0 0 


where Cy is the specific heat of the liquid at constant pressure, k is the thermal 
conductivity, and hp, denotes the latent heat of condensation. Writing this in 
terms of the boundary-layer variables, and letting AT = T,, — T.,, there isobtained 


w? 
9 Eh eteae a 
Fa, Gs hyp) = 2 O"(Ns) a) 
with Pr representing the Prandtl number (Py = v/a). Thus, a fixed value of 5 defines 
the parameter on the left-side of this equation. 


3. Approximate Solution 


The solution of the preceding equations can be obtained by numerically inte- 
grating equation (6) and matching this solution with equation (11) through use of 
the boundary conditions at 7 = 75. Equation (8) may then be solved through a 
separate numerical integration. Such a procedure would, however, be somewhat 
tedious in view of the number of parameters involved (0,/o, »,/», ys, Pr), and an 
approximate solution will instead be considered. 

It may readily be verified that series solutions of equations (6), (7), and (8) are 
of the form 


fn) = (>) [F i = 7 n° + | ; (13) 

" (é& Ge 
F(q) = Gy + Se! (n—y)® = 5 88 (Vt, (14) 
O(n) =1— 0+ 2S" SF ge... (15) 


where equation (14) follows from equations (10) and (11), Gj, = G@"(0), and # and o 
are constants. For all practical purposes, it will be necessary to evaluate equation 
(14) only in the neighborhood of ngs so as to satisfy the boundary conditions at 
7 = ng. The assumption will now be made that 75 and (75 — y) are sufficiently small 
such that the preceding equations may be rewritten as 


{(n) = (-) - - ume 
F(n) = Gy + 2! (9-9), (17) 
O(n) =1—o7n. (18) 


Physically, the neglect of higher-order terms in equation (13) is the same as 
neglecting acceleration effects within the film, since equation (16) gives 07u/dy =i); 
In turn, equation (18) neglects thermal convection within the film, and the resulting 


temperature profile is linear. 
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Equations (16), (17) and (18) contain five unknowns (B, v, o, G,, Gj), and four 
of these may be evaluated from equations (9b) giving o = 1/y5 and 


BS (2 (22 iy Gry (19a) 
0 y w 
W2 
|e (19) 


Gg — Sw (% (22) (2+)"] 2 19 
ao 1 — 5° ec 

w 2 ( y ) y 0 5 ( ) 

The remaining unknown, G/, has been tabulated as a function of G,, by EMMons 


and LeicH [4], and an abbreviated listing of their results is given in Table 1. 


Table 1. Values of G/) [4] 


_ 
ea nae | 0-369 | 0-406 | 0-483 | 0-563 | 0-645 ] 0-729 | 0-945 | 1-169 | 2-590 | 5-049 | 


w 


Now, from equations (16) and (18), equation (12) gives 


4/6, Ads 
rah pai 20 
Ns B (2 i) , ( ) 


Further, for most physical applications (0,/o <1) equation (19c) may be rewritten 


as G 
ey 
meh aa) ee Yao) 
G" x (-*) 8 


w 


and combining this with equations (19a) and (20) 


Gy z (4)28 { y \1/2 (& ( Cp AT 2/3 (21) 
(Gi)1i8 2 iG.) ae Pr i,,) 3 
Employing Table 1, the quantity G’ is completely determined for a given value of 
the parameter 


bre) elem 


Asymptotic expressions may also be obtained for G”. If the right-side of equation 
(21) is very small, then G,, ~ 0 and 


Gi, = 0:332. (22) 


On the other hand, if the right-side of equation (21) is large, then G, ~ 2G" /4] 
and equation (21) reduces to 


Coa em en leal ce |) a 


4, Results 


The results of usual practical importance are the shear stress and heat transfer 
at the plate surface. Consider first the local shear stress t, which is given by 


r=n(5") Ue Be) pro 


0 wot [0-2 | 04 | 0.6) [ose eo eee! 4... |)» Soui jae 
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From equations (16) and (19a), it is found that 
; Tv ? / u 
b 8 ) pte = Gy, (24) 
Oy Yoo 


where Re is the vapor Reynolds number defined as Re = Ux #/V 

The shear stress results are illustrated in Figure 2. For small values of the 
abscissa, a constant value is asymptotically approached given by equation (22). 
Equation (23) represents the asymptotic result for large values of the abscissa. 
One may note that since acceleration effects have been neglected the shear stress is 
constant throughout the film, and the results for t also apply at the liquid-vapor 
interface. 


' | — [ eee et 
} 
4 pecan Equation(23) | 
—-—fquation(22) | 
4 
A 
o 
== 
e 2 
JS 1 
: : 1 eer All ae SNe) 
0] 02 04 06 081 2 4 fey got YUU) 
1 
(x)? (2) (2247) fF. 
vw Oy Prey 
Figure 2 


Variation of local shear stress. 


With respect to heat transfer at the plate surface, the local heat flux per uni 
area q is given by : 


00 penn |ieaue, 
== = peel, nl ee) (2)\(0)) 7 
q= —h ToT) (5), = — (feo — Ta) |/ Fe 0°00) 


From equation (18), and writing the result in terms of a local Nusselt number 


a q * a 1 oe 
ieyt ee aE ae el 


Combining this with equations (19a) and (20), then 


se (a —(PEE GHy". es 


v 


A plot of the Nusselt number is shown in Figure 3. For small values of the 
abscissa, the asymptotic relation is found from equations (22) and (25) to be 


aie te Medea 2) me) 
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Upon combining equations (23) and (25) 


Nu (Zy" EOF (27) 
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which is the asymptotic expression for large values of the abscissa. 
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Variation of local Nusselt number. 


5. Discussion of Results 


It will be recalled that higher-order terms have been neglected in equation (13), 
and that this corresponds to negligible acceleration effects within the liquid film. 
The neglect of acceleration terms is therefore permissible if the condition ; 


2 


B 
“5, 3 <P 5 


is satisfied. From equation (20), this may be rewritten as 


In a similar manner, thermal convection within the film will be negligible only if 
higher-order terms are small in equation (15), and this gives 


CA 
ee SP 
Npg 


a 


A final requirement follows from equation (14) to be 


Therefore, the present analysis is applicable if the last two conditions are satisfied. 
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Admittedly, the present results neglect any ripple effects within the liquid film, 
and this could be an important factor in the actual physical process. Nevertheless, 


the analysis gives an insight into the physical phenomena involved, and illustrates 
the effect of the various parameters. 
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Zusammenfassung 


Die laminare Strémung von gesattigtem Dampf iiber eine gekiihlte ebene 
Flache wird untersucht fiir den Fall der Kondensation eines Films an der Ober- 
flache. Es wird angenommen, dass keine Gravitationskraft vorliegt und dass die 
Bewegung des Kondensats lediglich durch die Reibung am strémenden Dampf 
hervorgerufen wird. Ergebnisse fiir die Scherkraft und die Warmeiibertragung an 
der Oberflache werden mitgeteilt. Die wesentlichen Rechenergebnisse konnten in 
zwei Diagrammen zusammengefasst werden: Figur 2 gibt die zu erwartenden Schub- 
spannungen wieder, Figur 3 den Warmeiibergang. 


(Received: April 21, 1960). 


Signalfluss und Analogiedarstellung der Neutronenédkonomie 
beim Abbrand von Reaktorbrennstoffen 


Von WALTER HALG, Ziirich!). 


1. Allgemeine Bemerkungen 


Bei samtlichen Untersuchungen, welche die Gestehungskosten der Atomenergie 
zum Gegenstand haben, bildet die Frage nach der 6konomischen Verwendung der 
im Reaktor vorhandenen Neutronen eine ausschlaggebende Rolle. Jedes Neutron, 
welches aus dem Reaktorkern entweicht oder welches nicht zur Einleitung einer 
Kernspaltung oder zur Umwandlung von nichtspaltbarem Stoff in spaltbare Kerne 


- Verwendung findet, stellt ein Ausgabenposten in der Bilanz dar. Lewis?) hat ver- 


schiedentlich gezeigt, dass es zweckmAssig ist, die Neutronenbilanz der Kernbrenn- 
stoffe von derjenigen der iibrigen Reaktorbaustoffe einschliesslich Leckanteil 
separat zu betrachten. Man erreicht dadurch eine weitgehend universelle Behand- 


_ lungsweise der verschiedenen Reaktorkonzeptionen, weil sich der komplexe zeit- 
‘liche Ablauf von Neutronenproduktion und Neutronenkonsumation im Spaltstoff 


j gleichartig darstellt. 


i 
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1) ETH und Eidg. Institut fur Reaktorforschung, Wurenlingen. in . 
2) Vergleiche zum Beispiel W. B. Lewis, Lew cost fuelling without recycling, AECL, Nr. 3882, oder 


AECL Nr. 651. 
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Lewis definiert fiir den Einheitsneutronenfluss einen Neutroneniiberschuss o,(t) 
pro urspriinglich vorhandenem spaltbarem Kern Ny;.(0), fiir ein anfangliches Ge- 
misch aus U5 und U238 zum Beispiel 


o,(t) = No [2 N;, 06; € P — (2 N, 6; + (N an +so420+a2+)| . (4.1) 


fiss 


Dabei bedeuten N, die Konzentrationen der Atome eines der Figur 1.1 ent- 
sprechenden Spaltstoffzyklus?). 

Die Summation erstreckt sich iiber samtliche mit thermischen Neutronen spalt- 
baren Kerne, welche zum gewahlten Zyklus gehé6ren, so fiir den mit U** gestarteten 
Uran-Plutonium-Zyklus iiber i = 25, 49, 41 und beim Th*#?-U88-Zyklus tber 
j = 23, 25. 


Die Wirkungsquerschnitte*) sind dem speziellen Neutronenspektrum anzupas- 


sen; sie konnen im allgemeinen Fall irradiationsabhangige Gréssen sein. Der erste 
Summand entspricht der durch den Spaltstoff erzeugten Neutronenzahl, wahrend 
der zweite Term die im Spaltstoff absorbierten Neutronen beriicksichtigt. Zur Ver- 
einfachung wurde eine abgekiirzte, fiir sich selbst sprechende Schreibweise®) ein- 
gefiihrt. 


Die Definition (1.1) enthalt, wenn man das Neutronenspektrum des Reaktors — 


als gegeben betrachtet, ausser # keine Gréssen, welche durch die Konstruktion des 


Reaktors beeinflusst werden, sondern nur Materialkonstanten des Spaltstoffes. Ist — 
der Wert von o, positiv, dann muss in einem kritischen System dieser Uberschuss — 
durch spaltstoffremde Absorption (N o), und Oberflachenverluste (N o);,, kom- 


pensiert werden. 


Die spaltstofffremde Absorption enthalt auch den in den Regelorganen kompen- 
sierten Anteil des Neutroneniiberschusses. Weil nun aber fiir einen bestimmten — 
Reaktor die Kompensationsterme durch die Konstruktion gegeben sind, folgt aus — 
der Kenntnis von o, unmittelbar die maximal mégliche Betriebsdauer eines Reak- 


tors durch 


2. Differentialgleichungen der Atomkonzentrationen 


fisy(9) = (N 6), =F (N C) ae : ‘ (1.2) 


Zur Berechnung von o,(t) ist die Kenntnis der Konzentrationen der Atomkerne 


des Spaltstoffzyklus notwendig. Die entsprechenden Differentialgleichungen fiir 
den mit U* gestarteten Pu-Zyklus wie auch fiir den Th®3?-U233-Zyklus sind in 
Figur 2.1 und 2.2 zusammengestellt. 

Gio, OF, of und of, bedeuten dabei effektive Absorptionsquerschnitte, welche 
neben der thermischen Absorption auch die Resonanzprozesse der betreffenden 
Kerne beriicksichtigen. Solche Resonanzeinfange kénnen alle von verschiedenen 
spaltbaren Kernen erzeugte Neutronen erleiden, sofern sie auf die entsprechende 
Resonanzenergie abgebremst sind. Daher enthalt zum Beispiel N, og, neben einem 
thermischen Anteil Ngo, noch einen Zusatzterm N; 9,7; Bg;, wobei P,,; die 


Wahrscheinlichkeit dafiir beschreibt, dass ein schnelles Spaltneutron des i-ten 


; ) Die Indizes entsprechen der iiblichen Kennzeichnung, wobei anstelle der Doppelindizes 
einfache verwendet werden, sofern Missverstandnisse ausgeschlossen sind. 


) V enn keine besonde eB zeichnun i i i i i 
} s 7 g ange eben Ist, be h 
j V 3 : i r e geg zienen sich die Symbole auf die Ab 


°) y symbolisiert die Spaltprodukte. 
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Figur 2.1 
Differentialgleichungen fiir die Atomkonzentrationen des U-Pu-Zyklus. 
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Figur 2.2 
Differentialgleichungen fiir die Atomkonzentrationen des Th-U-Zyklus. 
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spaltbaren Kerns die Resonanzenergie des U88-Absorbers erreicht und dort auch 
eingefangen wird). 

Fiir natiirliche oder nur schwach angereicherte Uransysteme bleibt N,(t) nahezu 
konstant, so dass auch s; als konstant angesehen werden darf. Dasselbe gilt fiir 
Systeme mit grossem Anfangsgehalt an lpia eet : 

Anders verhalten sich Pu24° und U?34, deren Konzentrationen zeitabhangige 

Funktionen darstellen, so dass auch die entsprechenden § zeitlich variabel sind. 
In erster Naherung wird dieser Effekt dadurch beriicksichtigt, dass man $ propor- 
tional zu den entsprechenden Atomkonzentrationen setzt, wobei ein eventuell zeit- 
lich veranderliches Resonanzintegral, welches in die Ausdriicke fiir )(¢) resp. 
po4(t) eingeht, auch in den nunmehr nicht mehr linearen Differentialgleichungen 
noch unberiicksichtigt bleibt. 

Wir schreiben also zum Beispiel 


No of = No % + ps N, 0; 1; Box(t) & No % + Nol?) oF N, 6; 1; Bet). (2-4) 
- - 


v 


Die Wahrscheinlichkeit P*,(¢) ist auf ein Ny-Atom bezogen und bleibt, sofern man 
die Effekte im Resonanzintegral unberiicksichtigt lasst, zeitlich konstant, das heisst 


PF (2) = PI, = e Pe (1 — Py)* e- 08, (2.2) 


wobei der zeitliche Mittelwert (1 — p,)* ebenfalls pro Nj-Atom gegeben sei. 
Die gewahlte Darstellung erlaubt, fiir jeden Brutstoff’) (b) einen zeitlich varia- 
blen Spaltstoffregenerationsfaktor y,(t) zu bestimmen. 


_ pro Zeiteinheit gebildeter Spaltstoff N,(t) 09 (#) 
pro Zeiteinheit verbrauchter Spaltstoff 7 N;,(t) o; * 
i 


Yp(t) = 


Sind mehrere Brutstoffe vorhanden, so sei der Regenerationsfaktor des Zyklus 
gegeben durch 
v(t) = Dyo(t) - (2.4) 
b 


Die tblicherweise verwendeten, urspriinglichen Regenerationsfaktoren (initial con- 
version ratio, initial fuel regeneration factor) ergeben sich aus den Anfangswerten 
7»(0) resp. (0). 

Der allgemeine Ausdruck y(t) fiir den U-Pu-Zyklus ist in der Formelzusammen- 
stellung Figur 2.1 enthalten. 


Die Spaltprodukte sind durch drei Pseudoprodukte N, ,; p = 1, 2, 3 gemass 
% 


der Darstellungsweise nach Hurst’) beschrieben, wobei die sich rasch aufbauenden 


(2.3) 


| 


Produkte mit grossen Einfangsquerschnitten wie Xe und Sm separat beriicksichtigt 


worden sind), 


3, Signalflussdiagramme 


Losungen der im vorhergehenden Abschnitt diskutierten Differentialgleichungs- 
systeme konnen numerisch mit Digitalrechenmaschinen erhalten werden. Solche 


©) (l= P5)- 

7 ‘ i Pee 

) Der Brutstoff ist eine Substanz, welche durch Absorption von thermischen oder Resonanz- 
neutronen auf neues, spaltbares Material fiihrt, also zum Beispiel Th?#2, U288, [j234, py240, 


8) D. G. Hurst, Calculated Cross-Sections of Irradia ISS1 
: , k ted Fission P : . 
®) Vergleiche auch AECL Nr. 651. a 


v 
4 


Da. 


| 
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sind beispielsweise in AECL Nr. 651 wiedergegeben, wobei das dort verwendete 
Gleichungssystem fiir den U-Pu-Zyklus geringfiigige Unterschiede gegeniiber dem- 
jenigen der Figur 2.1 aufweist. 

Es ist aber auch méglich, die Differentialgleichungen naherungsweise zu lineari- 
sieren und das daraus durch Laplacetransformation erhaltene algebraische lineare 
Gleichungssystem entweder mit klassischen Methoden oder mittels eines Signal- 
flussdiagramms!°) aufzulésen. Die letztere Methode ist besonders dann niitzlich, 
wenn man die Auflésung nicht durchzufiithren gedenkt, sondern ein Kriterium zur 
Vereinfachung des Gleichungssystems sucht. Der Rechenaufwand wird namlich, 
besonders bei gekoppelten U-Th-Pu-Zyklen, ausserordentlich umfangreich, so dass 
es wiinschenswert ist, weniger bedeutsame Riickwirkungen von erst nach langer 
Bestrahlungsdauer auftretenden Atomkernen zu eliminieren, 


4. Auswertung durch Analogierechner 


Die Signalflussdarstellung legt eine Auflédsung der Gleichungssysteme mittels 
Analogierechengeraten nahe, obwohl selbstverstandlich fiir diese Auswertung die 
Erstellung eines solchen Diagramms keine Notwendigkeit bildet. Im allgemeinen 
wird eine Analogierechenmaschine mit Multiplikationseinheiten ausgeriistet sein, 
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Figur 3.1 
Signalfluss fiir Uran-Plutonium-Zyklus 


: 10) Uber die Erstellung des Signalflussdiagramms sowie die Regeln zu dessen Auflésung ver- 


- gleiche zum Beispiel S. J. Mason, Feedback-Theory; Further Properties of Signal Flow Graphes, 
“Proc. IRE 44 (7), 920 (1956). 
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Signalfluss fiir Thorium-Uran-Zyklus. 


so dass auch die im Signalfluss nicht zur Darstellung gelangenden nichtlinearen 
Terme mit beriicksichtigt werden kénnen. 

Figur 4.1 stellt das Analogieschaltbild des U-Pu-Zyklus dar, nach welchem mit 
der am Institut fiir Atomenergie in Kjeller, Norwegen, vorhandenen Maschine™) 
die nachfolgenden Resultate erhalten wurden. 

Obwohl die Differentialgleichungen (Figur 2.1) geringe Abweichungen gegen- 
iiber den von LEwis verwendeten aufweisen!”), ergibt eine Gegeniiberstellung 
(Figur 4.2) einiger in AECL Nr. 651 wiedergegebenen Resultate und jenen eines 
druckenden Digitalvoltmeters einen Hinweis auf die Genauigkeit unserer Analogie- 
rechnungen. 

In Figur 4.3 sind Registrierkurven fiir verschiedene Anreicherungen des ur- 
spriinglichen Spaltstoffes als Funktion der Irradiation (Neutronenfluss x Zeit) 
wiedergegeben. Dabei sei nochmals bemerkt, dass alle Gréssen mit Ausnahme des 
Regenerationsfaktors pro urspriingliches spaltbares Atom angegeben sind. Dies hat 
man bei Vergleichen, besonders bei der Bestimmung der maximalen Reaktor- 
betriebsdauer, zu beriicksichtigen. Die Terme fiir Fremdabsorption und Oberflachen- 
verluste sind ebenfalls in den entsprechenden Einheiten einzufiihren. 


a PACE 231 R (30 Integratoren, 45 Summationsverstarker, 100 Potentiometer, 4 Mehrfach- 
Multiplikatoren). 

Ich danke auch an dieser Stelle Herrn B. WINTHER sehr fiir die mir gebotene Gelegenheit, die 
Maschine zu benutzen und fiir seine tatkraftige Hilfe bei meinen Berechnungen. 


‘ a) Lewis verwendet einen mit der Irradiation veranderlichen Wirkungsquerschnitt fiir Plu- 
onium, 


> 
— 


Vol. XI, 1960 Kurze Mitteilungen — Brief Reports — Communications bréves J 


> 


[~] +Ng(O)o, 


oe t 
G7) 
LL @ e 
(Na); 0 e 
ma 


+ 
7 Ss Por 
aoe 
| _ és if 
A % 
fo jo i> 
6)— | N 
8D nama ae 
aaa i @3 
[=p @-) we Ser 
ad 
ba 2 fp Wols_ PEN, 0; pita oy a |b 
~ 2(NO)g (No), (No)y d 2(N0)p,, EN, 0,23 (NO)y42+6 (No)g +0 ty 
C | E F D 
0.5 V/mm 2V/mm 2V/mm 05 V/mm V/imn = = 1V/mm 
Figur 4.1 


Analogieschaltbild fiir Uran-Plutonium-Zyklus. 
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Figur 4.3 
Uran-Plutonium-Zyklus fiir verschiedenen U?-Gehalt. 
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Die Kurve des Regenerationsfaktors y durchlauft ein Minimum. Dieses tritt in- 
erster Naherung bei derjenigen Irradiation auf, wo die mit dem Wirkungsquer- 
schnitt multiplizierte Zerfallsrate von U*** gerade denselben Wert des entsprechen- ; 


den Produktes bei Pu®*® annimmt. 


Summary 


The Laplace transforms of the linearized equations for long term irradiation of | 
reactor fuel are arranged as a signal flow diagram which can either be used to 
solve the system or will make clear the correlation of important terms. 

Solutions for the Uranium-Plutonium cycle obtained with an analog computer 
are compared with Canadian results from digital evaluations. A general definition 
for the irradiation dependent fuel regeneration factor is introduced. ; 


(Eingegangen; 16, Mai bzw. 15. Juli 1960). 
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2nd ILMAC from 15th to 20th October 1962 


According to the great success of the first ILMAC, International Congress and 
Exhibition of Laboratory, Measurement and Automation Techniques in Chemistry, — 
which took place in November 1959 in Basle, the organizing committee decided 
unanimously to repeat that event. The second ILMAC will be accomplished from 
15th to 20th October 1962 again in the building of Swiss Industries Fair in Basle. 

M. TROTTMANN 
i 


Journal of Mathematical Analysis and Applications. The first number 
of the Journal of Mathematical Analysis and Applications has been released by 
Academic Press Inc., New York. Volume 1, Number 1, contains contributions by 
N. Levinson, A Nonlinear Volterra Equation Arising in the Theory of Superfluidity ; 
J.M. RicHarpson, The Hydrodynamical Equations of the One Component System 
Derived from Nonequilibrium Statistical Mechanics; H. Osporn, The Dirichlet Fune- 
tional. I; J. G. KemMENy, Ergodic Theorem for General Functions; R. P. Boas, Jr., 
Inequalities for Monotonic Series; A. ScHILD, On the Radiation Emitted by an Accele- 
vated Point Charge. 

The Editor of the Journal of Mathematical Analysis and Applications is Dr. dl 
BEeLLMAN, the RAND Corporation, 1700 Main Street, Santa Monica, California. 
Volume 1 will consist of four issues, to be released in 1960 at $16.—. Subscription 
orders should be sent to the Publishers, Academic Press Inc., New York 3, or 
London, S. W. 1, England. L. HENtEr 
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Demmndachst erscheint: 
LINDER 


Statistische Methoden fiir Naturwissenschafter, 
Mediziner und Ingenieure 


Von Prof. Dr. ARTHUR LINDER, Professor an der Universitat Genf und an der Eidgendéssi- 
schen Technischen Hochschule in Ziirich. 484 Seiten mit 58 Figuren und 73 Beispielen. 
Preis gebunden Fr./DM 54,— 

Dritie, umgearbeitete und stark erweiterte Auflage 
Mathematische Reihe, Band 3. Sammlung «Lehrbiicher und Monographien aus dem Gebiete 
der exakten Wissenschaften», 


Angesichts der unvermindert anhaltenden Nachfrage hat sich der Verfasser entschlossen, 
seine bekannte Einfiihrung in die Methoden der mathematischen Statistik stark zu erwei- 
tern. Dabei wurde gleichzeitig der den Anwendungen gewidmete Teil véllig umgearbeitet. 
Neu aufgenommen wurden unter anderem: Verschiedene Anwendungen von Chiquadrat, die 
Streuungszerlegung bei ungleichen Klassenhaufigkeiten, die Bestimmung von Streuungs- 
komponenten, die nichtlineare Regression, die «analysis of covariance», das Schatzen von 
Parametern und einige Transformationen yon Prozentzahlen (Arc sin, Probit, Logit, Loglog). 


Vom gleichen Autor erschien in unserem Verlag: 


Planen und Auswetten von Versuchen 


Eine Einfiihrung fir Naturwissenschafter, Mediziner und Ingenieure. 182 Seiten mit 

9 Figuren. Ganzleinen Fr./DM 21.—. Zweite Auflage 1959. (Reihe der experimentellen 

Biologie, Band 13 —- Sammlung «Lehrbiicher und Monographien aus dem Gebiete der 
exakten Wissenschaften ».) 


Diese Einfiihrung richtet sich an Naturwissenschafter, Mediziner und Ingenieure; sie setzt 
keine Kenntnisse in mathematischer Statistik voraus. Der Leser wird angeleitet, Versuche 
richtig zu planen und einwandfrei auszuwerten. Die griindlich durchgearbeiteten Anwen- 
dungsbeispiele aus der biologischen, medizinischen, industriellen und landwirtschaftlichen 
Forschung bilden einen wichtigen Bestandteil des Buches. Mit besonderer Sorgfalt wurden 
die grundlegenden Gedankengange herausgearbeitet. Wegleitend fiir die Auswahl und 
Anordnung des Stoffes waren die Erfahrungen, welche der Verfasser aus Vorlesungen 
schépfte, die er seit 1938 an der Universitat Bern, an der Eidg. Technischen Hochschule 
in Ziirich und an der Universitat Genf hielt. 4 


1961 erscheint: 


Handliche Sammlung 


mathematisch-statistischer Tafeln 
40 Seiten, Fr./DM 4.50 


Die vorliegende Sammlung statistischer Tafeln soll die umfangreicheren, bisher erschiene- 

nen nicht ersetzen: sie stellt vielmehr einen handlichen Auszug der am meisten verwen- 

deten Tafeln dar, die dazu bestimmt sind, stets griffbereit auf dem Arbeitstisch zu legen. 

Ausser den Tafeln von ¢, Y? und F enthalt diese Sammlung Tafeln fiir die Transforma- 

tionen von Prozentzahlen mittels Arcussinus, Probits, Logits und Loglog, sowie Tafeln 

fiir zufallig angeordnete Zahlen, die besonders beim Planen von Versuchen und Stich- 
probenerhebungen zu beniitzen sind. 


Zu bexichen durch Ihre Buchhandlung — Obtainable from your bookseller — 
Commandes 4 votre libraire 
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